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PARAMETRIC  INVESTIGATION  OF  RADOME 
ANALYSIS  METHODS:  SALIENT  RESULTS 

I .  Introduction 

The  final  technical  report  on  this  three-year  radome  research  program 
consists  of  four  volumes.  This  report,  Volume  I,  presents  the  talient  results 
conclusions,  and  recommendations  of  this  research  whose  objective  was  to  deve¬ 
lop  a  general  theory  of  radome  analysis  and  to  determine  the  accuracies  and 
ranges  of  validity  of  three  particular  computer-aided  radome  analysis  methods. 
The  results  are  presented  in  appendices  as  copies  of  papers  prepared  for  publi 
cation  as  summarized  in  the  main  body  of  this  report. 

Volume  II  documents  the  analytical  method  and  Fortran  computer  code 
used  to  analyze  the  antenna/radome  combinations  using  a  fast  receiving  for¬ 
mulation  based  on  Lorentz  reciprocity  and  geometrical  optics.  Volume  III  docu 
ments  the  analytical  method  and  additional  Fortran  software  required  for  ra¬ 
dome  analysis  based  on  the  Huygens-Fresnel  principle  or  surface  integration. 
Volume  IV  documents  the  pattern  and  bore sight  error  measurements  made  on  eight 
combinations  of  three  monopulse  antennas  and  five  tangent  ogive  radomes  at 
35  GHz. 

Measured  data  is  used  as  true  data  in  assessing  the  accuracies  of  the 
computer  codes.  It  is  expected  that  the  measured  data  obtained  will  be  used 
in  the  future  for  similar  purposes.  It  is  also  expected  that  the  documented 
computer  codes  will  serve  as  part  of  a  technology  base  for  use  by  researchers 
and  practitioners  in  the  radome  technical  community. 

II.  Results 

Appendix  A  presents  a  paper  which  describes  the  theory  of  radome  analy¬ 
sis  developed  during  this  investigation.  The  theory  is  simply  an  application 
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of  the  Huygens- Fresnel  principle  and  Lorentz  reciprocity  to  the  radome  pro¬ 
blem.  All  methods  of  radome  analysis  appearing  in  the  literature  are  embodied 
in  this  general  theory.  It  provides  the  correct  framework  in  which  to  think 
about  radome  analysis  and  to  make  objective  comparisons  between  the  various 
methods  of  analysis. 

Appendix  B  presents  a  paper  describing  some  antenna  synthesis  work  that 
was  undertaken  during  this  research.  A  procedure  is  described  whereby  the 
aperture  fields  of  the  four-horn  monopulse  antenna  without  radome  can  be 
synthesized  from  measured,  amplitude -only,  far-field,  principal  plane  patterns. 
A  priori  information  about  the  geometry  and  excitation  of  the  actual  antenna 
is  used  to  find  a  credible  solution  to  an  inverse  source  problem  for  which 
there  would  otherwise  be  no  unique  solution.  In  the  absence  of  the  radome, 
each  analysis  method  accurately  predicts  the  measured  patterns. 

Appendix  C  presents  comparisons  of  measured  and  computed  boresight 
«rrors  and  radome  losses  for  four  of  the  eight  antenna/radome  combinations 
considered.  The  antennas  and  radomes  used  are  briefly  described.  (More  de¬ 
tailed  descriptions  are  presented  in  Volume  IV.)  Qualitative  conclusions 
about  the  accuracies  of  the  methods  can  be  made  from  the  data  presented;  in 
addition,  recommendations  for  future  work  are  brought  to  light. 

III.  Conclusions  and  Recommendations 

From  the  data  presented  in  Appendix  C,  it  is  concluded  that  none  of 
the  computer  aided  radome  analysis  methods  investigated  consistently  and 
accurately  predicted  the  boresight  errors  and  losses  of  the  various  antenna/ 
radome  combinations  used.  For  moderate  size  antennas  and  radomes  (as  defined 
in  Appendix  C) ,  the  fast  receiving  method  predicts  boresight  errors  reasonably 
accurately  and  is  perhaps  the  fastest  method  available  anywhere;  however,  for 
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small  radomes  and  large  radomes,  the  accuracy  of  the  error  predictions  is 
not  as  good. 

The  surface  integration  method  was  not  completely  assessed  for  pre¬ 
dictive  accuracy  because  of  unresolved  problems  with  the  code.  A  large  part 
of  tne  problem  with  the  development  of  the  code  was  the  relatively  long  exe¬ 
cution  times  required  on  the  Cyber  70  system  used.  It  is  clear  that  because 
of  the  long  execution  times,  use  of  the  surface  integration  method  will  be 
restricted  to  the  analysis  of  small  radomes. 

It  is  recommended  that  some  "fine  tuning"  be  done  on  the  fast  receiving 
method  to  improve  its  predictive  accuracy  over  a  larger  range  of  antenna/ra- 
dome  parameters.  Specifically,  the  ray  tracing  procedure  should  be  modified 
to  account  for  refractive  effects  and  propagation  of  the  fields  along  the  rays. 
The  expected  improvement  in  accuracy,  coupled  with  the  inherent  practicality 
of  this  code,  make  this  recommended  endeavor  worthwhile. 

It  is  also  recommended  that  work  be  continued  on  the  development  of 
the  surface  integration  method  of  analysis  so  that  its  predictive  accuracy 
can  be  clearly  established.  A  dedicated,  small  computer  system  (32-bit  word 
length)  is  recommended  for  this  work  to  avoid  the  anticipated  high  costs 
associated  with  very  heavy  use  of  a  time-share  system  such  as  the  Cyber  70 
at  Georgia  Institue  of  Technology. 

It  is  recommended  that  exact  solutions  of  carefully  selected  antenna/ 
radome  configurations  be  obtained  for  use  as  true  data  in  the  assessment  of 
the  accuracies  of  various  computer-aided  radome  analysis  methods.  Configura¬ 
tions  used  must  conform  to  those  expected  in  the  applications.  The  accuracies 
of  the  solutions  and  the  solutions  themselves  must  be  impeccable.  Until  such 
time  as  this  recommendation  is  carried  out,  measured  data  must  continue  to  be 
used  as  true  data. 
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As  a  parallel  approach  to  exact  solutions,  it  is  recommended  that 
numerical  solutions  be  obtained  on  practical  configurations  of  interest 
using  modern  numerical  methods  such  as  method  of  moments  and  integral  equa¬ 
tion  formulations.  The  use  of  specialized  computing  structures  for  dedicated 
application  to  such  elec tr^maanc tic  radiation  and  scattering  problems  should 
be  investigated. 

Finally,  it  is  recommended  that  experimental  techniques  be  developed 
to  help  isolate  deficiencies  in  the  methods  of  analysis.  The  usual  measui.  ■ 
ments  of  patterns  and  boresight  errors  simply  do  not  provide  the  necessary 
information  to  pinpoint  invalid  assumptions  and  poor  approximations  in  any 
analysis  procedure.  New  experimental  methods  which  make  use  of  automated 
measurements  and  near-field/far-field  transformations  need  to  be  developed. 


4 


APPENDIX  A 


"Theory  of  Radome  Analysis" 


Submitted  for  Review  for  Publication 


in 


IEEE  Transactions  of  the  Antennas  and  Propagation  Society 


November  1980 


5 


THEORY  OF  RADOME  ANALYSIS 


G.  K.  Huddleston 
School  of  Electrical  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 

Abstract 

A  basic  theory  of  radome  analysis  is  presented  based  on  the 
reciprocity  theorem  and  the  Huygens-Fresnel  principle.  Receiving  and 
transmitting  formulations  are  developed.  Techniques  of  analysis  are 
presented  to  distinguish  the  salient  features  of  radome  analysis  and  to 
show  the  relationships  between  the  theory  and  some  existing  methods  of 
analysis.  The  equivalence  between  the  receiving  and  transmitting  cases 
is  also  established,  both  in  general  and  explicitly  for  the  far-field 


case . 
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I .  INTRODUCTION 

Radome  analysis  is  the  application  of  electromagnetic  theory  to 
determine  the  effects  of  protective  dielectric  structures  on  the  elec¬ 
trical  characteristics  of  antennas  enclosed  by  them.  Numerous  methods 
of  analysis  have  been  developed  for  the  prediction  of  radome  effects; 
however,  no  unified  theory  has  been  advanced  to  provide  a  common  basis 
of  understanding  of  the  various  approaches  taken  and  the  approximations 
used  therein.  Furthermore,  no  comprehensive  measurements  have  been  re¬ 
ported  which  provide  true  data  in  determining  the  accuracies  of  the 
various  methods  when  parameters  such  as  antenna  size,  wavelength,  ra¬ 
dome  size  and  shape  are  considered  over  the  ranges  that  they  are  likely 
to  assume  in  the  applications. 

This  paper  presents  a  theory  of  radome  analysis  which  embodies 
known  methods  of  analysis  that  have  appeared  in  the  literature.  The 
theory  is  based  on  the  reciprocity  theorem  [1]  and  the  Huygens-Fresnel 
principle  for  electromagnetic  fields  [2,3]  both  of  which  are  derivable 
from  Maxwell's  equations  via  the  divergence  theorem  and  .re  vector 
Green's  theorem.  The  reciprocity  theorem  serves  as  the  basis  for  the 
receiving  formulation  of  radome  analysis,  and  the  Huygens-Fresnel  prin¬ 
ciple  is  the  basis  for  the  transmitting  formulation  as  seen  in  what 
follows.  The  equivalences  of  the  two  formulations  are  also  established. 
Some  approximations  and  techniques  used  in  implementing  the  analyses  are 
presented  along  with  discussions  which  clarify  the  justifications  for 
the  approaches  used  in  practice. 
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The  development  and  presentation  of  this  theory  of  radome  analy¬ 
sis  is  motivated  by  an  on-going  parametric  investigation  of  radome  analy¬ 
sis  methods  [4].  Measured  pattern  data  and  boresight  err cr  data  for 
three  antenna  sizes  combined  in  fifteen  combinations  with  five  radome s 
have  been  obtained  to  serve  as  true  data  in  determining  the  accuracies 
and  ranges  of  validities  of  three  common  computer-aided  analyses.  Great 
care  has  been  taken  to  ensure  accurate  modeling  of  the  antennas  such  that, 
in  the  absence  of  a  radome,  each  computer  code  accurately  predicts  the 
measured  radiation  patterns  [5 ] .  The  theory  provides  the  framework  in 
which  the  methods  can  be  compared,  especially  in  regard  to  the  approxi¬ 
mations  used,  the  validity  of  various  assumptions  made,  and  the  equiva¬ 
lence  of  different  computational  procedures  which  yield,  or  should  yield, 
the  same  result.  The  theory  also  provides  the  basis  for  developing  test 
cases  to  verify  the  computer  codes. 

Development  of  radome  analysis  methods  have  paralleled  the  develop¬ 
ment  of  airborne  radar  systems  whose  complexities  have  increased  over  the 
years.  Silver  [6]  illustrates  the  geometrical  optics  approaches  taken  ur 
to  1949  and  which  were  developed  during  the  previous  war  years.  Kilcoyne 
[7]  presented  a  two-dimensional  ray  tracing  method  for  analyzing  radomes 
which  utilized  the  digital  computer,  and  is  an  extension  of  work  done 
earlier  by  Snow  [8]  and  by  Pressel  [9],  A  more  rigorous  method  of  analy¬ 
sis  was  introduced  in  the  same  year  by  Van  Doeren  [10]  using  an  integral 
equation  to  compute  fields  inside  the  radome  caused  by  an  incident  plane 
wave.  Tricoles  [11]  formulated  a  three-dimensional  method  of  radome 
analysis  based  on  Shelkunoff's  induction  and  equivalence  theorems.  Tavis 
[121  describes  a  three-dimensional  ray  tracing  technique  to  find  the 
fields  on  an  equivalent  aperture  external  to  an  axially  symmetric  radome. 


Paris  [13]  describes  a  three-dimensional  radome  analysis  wherein  the  tan¬ 
gential  fields  on  the  outside  surface  of  the  radome  due  to  the  horn  an¬ 
tenna  radiating  inside  the  radome  are  found.  Wu  and  Rudduck  [14]  describe 
a  three-dimensional  method  which  uses  the  plane  wave  spectrum  representa¬ 
tion  to  characterize  the  antenna.  Joy  and  Huddleston  [15]  describe  a 
computationally  fast,  three-dimensional  radome  analysis  which  utilizes 
the  plane  wave  spectrum  (PWS)  representation  and  exploits  the  Fast  Fourier 
Transform  (FFT)  to  speed  up  the  computer  calculations .  Chesnut  [16]  has 
combined  the  program  of  Wu  and  Rudduck  with  the  work  of  Paris  to  form  a 
three-dimensional  radome  analysis  method.  Huddleston  [17]  has  recently 
developed  a  three-dimensional  radome  analysis  method  which  uses  a  general 
formulation  based  on  the  Loren tz  reciprocity  theorem.  Siwiak,  et  al . , 

[18]  have  recently  applied  the  reaction  theorem  to  the  analysis  of  a  tan¬ 
gent  ogive  radome  at  X-band  frequencies  to  determine  boresight  error. 
Hayward,  et  al.,  [19]  have  compared  the  accuracies  of  two  methods  of 
analysis  in  the  cases  of  large  and  small  radomes  to  show  that  ray  tracing 
does  not  accurately  predict  wavefront  distortion  in  the  case  of  small 
radomes . 


II.  RECEIVING  FORMULATION 

The  reciprocity  theorem  is  a  starting  point  for  the  formulation 
of  a  unified  theory  of  radome  analysis.  The  general  reciprocity  theorem 
[1]  states  that 

<p  (E  x  H  -  K  xH  }*n  da  = 

Js  -a  -b  -b  --a 

(1) 


(E,  •  J  -  E  •  J.  -  H.  •  J10  +  H  •  J™)  *dV 

— b— a— a-b— b— a— a— b 
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X-- 


where  the  surface  S  encloses  the  volume  V  containing  two  sets  of  electric 

and  magnetic  sources  (J  ,  J*")  ,  (J,  ,  jf1)  which  give  rise  to  electromagnetic 

—a  — a  — b  b 

fields  (E  ,  H  ) ,  (E,  ,  H.)  ,  respectively.  (Time  variations  of  the  for r 
—a  —a  — b  — b 

e^wt  are  understood  and  suppressed)  . 

Figure  1  illustrates  a  typical  radome  analysis  situation  where  it 
is  desired  to  determine  the  response  V  of  Antenna  "a"  to  fields  produced 

cl 

by  Antenna  "b".  The  closed  surface  S  in  the  reciprocity  theorem  is  chosen 
to  be  Sq  +  so  that  the  volume  V  is  the  source-free  region  lying  between 
Antenna  "a"  and  the  inner  surface  of  the  dielectric  enclosure  depicted  in 
Figure  1.  In  this  case,  the  right  side  of  Equation  (1)  equates  to  2ero 
(Lorentz  reciprocity  theorem) ,  and  the  left  side  can  be  separated  into 
two  surface  integrals  as 


(E  x  H,  -  E,  x  H  )  *n  da  =  (E  xH  -  E.  x  H  )‘n  da  (2) 

-a  -b  ^b  — a  )_  -a  -b  -a 

so  S1 


where  n  is  the  unit  normal  to  the  surface  and  is  directed  positively  out 

of  the  volume  V.  The  fields  E  ,  H  are  those  produced  on  S,  and  s„  by 

—a  —a  10 

Antenna  "a"  when  it  is  activated  (transmitting) ;  E^,  are  the  fields 
produced  on  and  when  Antenna  "b”  is  activated  (and  Antenna  "a"  is 
passive  or  receiving) . 

I 

Let  surface  be  divided  into  two  parts:  consisting  of  that 

portion  of  SQ  across  the  waveguiding  structure  connecting  the  generator 

»» 

of  Antenna  "a"  to  the  radiating  structure;  consisting  of  the  remainder 
of  Sq  and  which  coincides  with  the  conducting  surface  of  the  antemia  as 

II 

indicated  in  the  figure.  The  integral  of  the  fields  over  is  identi¬ 
cally  zero  [20]  so  that  the  left  side  of  Equation  (2)  reduces  to  the 

I 

integration  over  S  .  Now,  in  the  many  practical  cases  of  interest, 

I 

there  exists  a  dominant  mode  of  propagation  over  (e.g.,  TE^Q 
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A 


FIGURE  1.  GEOMETRY  FOR  RECEIVING  FORMULATION 
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mode  in  rectangular  waveguide)  so  that  currents  and  voltages  can  be 
defined  at  this  terminal  plane  so  that  there  results 

' 

Vb  +  vb  Xa  =  -  *  4  -  4  X  V  *n  da  (3) 

S0 

The  currents  and  voltages  are  defined  for  the  two  cases  of  in¬ 
terest  in  Figure  2: 

I 

Va,  Ia  -  Voltage  and  current  produced  at  by  generator  "a" 

when  generator  "b"  is  passive; 

I 

Vfa,  1^  -  Voltage  and  current  produced  at  by  generator  "b" 

when  generator  "a"  is  passive. 

Linear  impedance  relationships  can  also  be  defined  as  indicated  in  Figure 
2  such  that  (using  Y  =  1/2) 


V  I  +  V  I 
a  b  b  a 


v  (V.  Y  )  +  v.  (V  Y.) 

cl  £)  3  D  cl  1 


Equating  Equation  (4)  to  Equation  (2)  and  solving  for  the  received  vol¬ 
tage  V^  results  in 


V  =  _ ^c_ 

b  V  (Y  +Y  ) 
a  a  1 


WV 


x  H.  -  E,  x  H  )  *n  da 
-b  -b  — a 


where  V^,^  has  been  defined  for  convenience.  This  is  one  desired  result. 

In  most  analysis  situations,  Antenna  "a"  and  the  dielectric  enclo¬ 
sure  comprise  the  antenna/radome  combination.  The  separation  of  this  com¬ 
bination  and  Antenna  "b"  can  be  made  very  large  so  that  the  fields  of  "b" 
are  those  of  a  transverse  electromagnetic  plane  wave  incident  on  the  ra- 
dome.  In  radar  applications,  the  distant  target  is  illuminated  either 
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(a)  Source  "a"  active, 


Source  "b"  passive 


(b) 


FIGURE  2. 


Source  "a"  passive, 


Source  "b"  active 


DEFINITIONS  OF  CIRCUIT  PARAMETERS  OF  THE  ANTENNAS. 
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by  Antenna  "a"  or  another  antenna  not  shown  in  Figure  1.  The  reflec¬ 
tion  from  the  target  comprises  Antenna  "b"  in  the  analysis  and,  for 
large  separation,  arrives,  again,  as  an  incident  plane  wave  on  the  ra- 
don<_ .  Note,  however,  that  the  receiving  formulation  embodied  in  Equa¬ 
tion  (r')  is  valid  for  either  the  near-field  (small  separation)  or  far- 
fi«=ld  (large)  case.  Also,  superposition  applies  so  that  multiple  tar- 
uets  or  sources  can  be  accommodated  as  linear  combinations  comprisino 

the  fields  E  .  H,  in  Equation  (5). 

— t)  — b 

Equation  (S)  is  exact;  however,  in  the  evaluation  of  the  inte¬ 
gral,  approximations  to  the  actual  fields  on  may  be  introduced  so 
that  inaccuracies  in  the  computation  of  the  received  voltage  may  result. 
To  see  the  effects  of  certain  approximations,  it  is  advantageous  to  make 
the  following  definitions  of  primary  and  scattered  fields: 


E  =  E1  +  E  +  E’  (6) 

—a  — ao  —as  —as 


H 

—a 


+  H 
—as 


+  H 
—as 


(71 


E  , 

— ao 


E 

—as 


H 

—as 


Fields  on  due  to  Antenna  "a"  which  contain  no 
scattered  field  components;  i.e.,  the  fields  that 
would  exist  cn  S  in  the  absence  of  the  dielectric. 
First  order  scattered  fields  on  due  to  Antenna  "a” 
i.e.,  fields  reflected  only  once  from  material  media 
outside  and  from  Antenna  "b". 

Higher  order  scattered  fields  on  due  to  Antenna 
■’a";  i.e.,  components  caused  by  multiple  scattering 
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between  Antenna  "a"  and  surrounding  media,  including 
Antenna  "b". 

Similar  definitions  can  be  made  for  Antenna  "b"  with  respect  to 
the  surface  as  follows: 


E.  =  El  +  E,_  +  E, 

—b  —bo  — bs  — bs 


(8) 


h  = 


+  Ik 

— bs 


(9) 


where  E^  ,  are  the  fields  on  S,  due  to  Antenna  "b"  and  which  contain 

-bo  -bo  1 

no  scattered  components  due  to  material  media  inside  S^,  including  Antenna 

"a".  The  fields  E  ,  H,  are  first  order  scattered  fields  on  S,  due  to 
-bs  -bs  1 

I  I 

a  single  scattering  from  Antenna  "a".  The  fields  E^,  are  the  higher 
order  scattered  fields.  When  these  definitions  are  substituted  into  Equa¬ 
tion  (5)  and  the  vector  cross  products  evaluated,  there  results 


V„„  =  l_  (E1  X  h5  -  e5  x  H1  )  *n  da 
REC  JS^  — ao  —bo  —bo  — ao 


_  (E1  x  H.  -  E,  x  H1  )  *n  da 
— ao  -bs  — bs  — ao 


+|  (E  x  -  e5  x  H  )*n  da 
J  ~as  —bo  —bo  —as 

1 


(E  x  H  -  E,  x  H  )*n  da 
-as  -bs  -bs  -as 


+  Terms  involving  higher  order  scattered  fields  (10) 


The  first  integral  on  the  right  side  of  Equation  (10)  provides 
the  primary  contribution  to  the  received  voltage  for  those  cases  of  in¬ 
terest  where  the  dielectric  is  somewhat  transparent  to  the  incident 

fields  E.1  ,  H,1  of  Antenna  "b".  The  second  and  third  integrals  are  each 
—bo  —bo 

identically  zero  by  virtue  of  the  definitions  of  the  scattered  fields  and 
the  general  reciprocity  theorem  as  applied  to  the  source-free  volume  V 
consisting  of  the  region  outside  S^;  i.e.,  since  there  are  no  sources  of 
the  defined  fields  in  this  region,  the  right  side  of  Equation  (1)  and, 
hence,  the  second  integral  of  Equation  (10) ,  are  both  zero.  A  similar 
argument  holds  for  the  third  integral  when  the  general  reciprocity  theo¬ 
rem  is  applied  to  the  source-free  region  V  enclosed  by  . 

The  contribution  of  the  fourth  integral  in  Equation  (10)  depends 
on  the  scattering  properties  of  Antenna  "a"  as  well  as  on  the  reflec¬ 
tive  properties  of  the  dielectric  enclosure  with  respect  to  the  incident 

fields  E1  ,  H1  .  It  is  difficult  to  assess  the  significance  of  these 
— ao  — ao 

contributions  and  those  of  the  integrals  involving  higher  order  scattered 
fields.  At  any  rate,  it  is  impractical  in  current  computer-aided  analy¬ 
ses  to  include  scattering  from  the  antenna  inside  the  enclosure;  hence, 
only  the  first  term  in  Equation  (10)  is  usually  retained.  The  objective 
of  practical  radome  analysis  then  centers  on  the  determination  of  the 
fields  in  the  integral 


(E1  x  -  E^  x  H1  ) *n  da  ( 1 ) ) 

_  — ao  —bo  -^o  — ao 

S1 

and  on  the  evaluation  of  the  integral  itself. 

Surfaces  other  than  the  inside  surface  of  the  dielectric  of 

Figure  1  can  also  be  selected  for  the  evaluation  of  the  received  voltage 
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I 


via  Equation  (10)  or  (11) .  For  example,  choose  SQ  as  before  and  choose 

to  also  coincide  with  the  outer  surface  of  Antenna  "a",  including  the 

radiating  aperture  portion  S  .  Apply  Lorentz  reciprocity  to  the  source- 

ap 

f 

free  volume  contained  between  the  antenna  aperture  and  to  yield 


ss  <e1  x  h5  -  e5  x  H1  )  *n  da  (12) 

REC  Js  “ ao  -bo  —bo  — ao 

ap 

The  fields  E1  ,  H1  are  the  aperture  fields  of  Antenna  "a"  when  it  is 
-ao  -ao 

transmitting  in  the  absence  of  the  dielectric  and  are  often  known  or  spe¬ 
cified  for  analysis  purposes.  The  real  difficulty  arises  in  the  determi¬ 
nation  of  the  fields  ,  h5  on  the  chosen  surface  S,»  And  it  is  the 

—bo  —bo  1 

choice  of  the  surface  and  the  approximations  used  to  find  the  fields 
on  it  which  distinguish  the  various  methods  of  analysis  based  on  the  re¬ 
ceiving  formulation.  More  will  be  said  concerning  the  determination  of 
these  fields  in  the  next  section. 

Another  choice  of  the  surface  of  integration  consists  of  the 
outer  surface  of  the  dielectric  enclosure  of  Figure  1.  First  apply 
Equation  (1)  to  the  source-free  region  consisting  of  the  dielectric  itself 
as  enclosed  by  its  inner  surface  and  its  outer  surface  S^;  i.e., 

S  *  +  S^.  Since  the  right  side  of  Equation  (1)  is  zero,  and  since  the 

positive  direction  of  n  is  out  of  V,  it  is  seen  from  Equation  (5)  that 


REC 


(E 


x  H  )'n  da 


(13) 


where  n  is  the  unit  outward  normal  to  the  outer  surface  of  the  dielectric. 
It  will  be  shown  in  a  later  section  that  Equation  (13)  is  equivalent  to 
the  transmitting  formulation  to  be  discussed  next. 
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III.  TRANSMITTING  FORMULATION 


A  second  generalized  approach  to  radome  analysis  uses  a  transmit¬ 
ting  formulation  which  does  not  explicitly  consider  the  fields  produced 
by  the  incident  fields  of  a  source  outside  the  dielectric  structure. 
Instead,  the  equivalence  and  uniqueness  theorems  of  electromagnetics  [21] 
are  invoked  to  establish  the  result  that  the  fields  radiated  by  Antenna 
"a"  into  the  unbounded,  homogeneous  region  outside  the  dielectric  en¬ 
closure  can  be  found  from  knowledge  of  the  tangential  electric  and  mag¬ 
netic  fields  on  the  outer  surface  of  the  dielectric.  Moreover,  there 
exists  an  integral  formula  for  the  actual  computation  of  these  fields. 

After  Stratton  and  Chu  [2],  Silver  [3]  derives  the  general  solu¬ 
tion  for  the  time-harmonic  electromagnetic  fields  in  a  homogeneous  medium 
which  arise  from  a  prescribed  set  of  sources,  including  magnetic  (equiva¬ 
lent)  charges  and  currents.  The  derivation  is  based  on  Maxwell's  equa¬ 
tions  (including  the  equations  of  continuity  for  charges  and  currents) , 
a  vector  Green's  theorem,  and  the  free-space  Green's  function  i|>  =  e 
The  results,  Equations  (3.108)  and  (3.109)  of  Silver,  are  called  the 
Huygens-Fresnel  principle  (alsr>  called  Kirchhof f-Huygens  principle)  and 
are  repeated  here  for  convenience: 


E (x ' , y ' ,z ' ) 


(~jwyi))(n  x  E)  +  (n  x  E)  xV^  +  (n*E)Vi/*]  dS 
S 


(14) 


H (x ' ,y ' , z ' ) 


T 

4tr 


[ j <ue (n  x  EH  +  (n  x  H)  xV\p  +  (n*H)V^]  dS 
S 


(15) 


where  n  is  the  unit  normal  to  S  which  points  into  the  source-free  medium 
as  indicated  in  Figure  3. 
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In  the  above  equations,  primed  variables  are  the  coordinates  of 


the  point  at  which  the  field  is  to  be  computed,  and  unprimed  coordinates 
are  used  to  designate  the  source  coordinates  lying  on  the  surface  S.  The 
distance  r  is  measured  from  the  field  point  (x',y’,z’)  to  the  source  point 
(x,y,z)  so  that 


r  =  Ax-x')*  +  (y-y ' )  /  +  (z-z')z 


(16) 


and 


Vtp  =  -  ( jk  +—■)<()  r  =  (jk  +  ^-)  ip  r  ^  (17) 

where  r  =  -r^  is  a  unit  vector  directed  along  r  from  the  field  point  to 
the  source  point.  The  variable  T  is  adopted  from  the  notation  of  Poggio 
and  Miller  [22]  so  that  the  fields  at  points  on  S  can  also  be  included; 
it  is  given  by 


T 


1 

l-Q/4ff 


1  if  n=0,  (x',y’,z’)  not  on  S 

2  if  (1=2tt,  (x',y',z')  on  S 


(18) 


The  terms  (n  x  E)  ,  (n  x  H)  in  Equations  (14)  and  (15)  can  be  con¬ 
sidered  as  equivalent  magnetic  and  electric  surface  current  densities 


K  =  n  x  H 


(19) 


K  =  -  (n  x  E) 

TD  — 


(20) 
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with  corresponding  equivalent  surface  charge  densities 


rj  =  (n  .  E) 


(21) 


h  =  (n  •  H) 
m  — 


(22) 


These  equivalent  currents  and  charges  are  the  sources  of  the  fields  E,  H 
in  the  region  outside  the  dielectric  enclosure.  Using  the  equation  of 
continuity  for  surface  currents  and  charges  [23] . 


V 


S 


K  + 


la 

at 


o 


(23) 


the  terms  in  Equations  (14)  and  (15)  involving  normal  field  components 
may  be  rewritten  for  the  time-harmonic  case  as 


n  •  E  =  V  •  (n  x  H)  (24) 

—  U )E  S  — 

n  •  H  =  ^-  Vg  •  (E  x  n)  (25) 

where  V  •  is  the  "surface  divergence".  Thus,  the  fields  E,  H  in  V  are 
expressible  entirely  in  terms  of  the  tangential  fields  on  S  as  is  in  keep¬ 
ing  with  the  uniqueness  theorem.  (Note,  however,  that  the  surface  diver¬ 
gence  involves  first  order  partial  derivatives  whose  numerical  computa¬ 
tion  may  not  be  desirable.  ) 

Implicit  in  the  application  of  equivalent  surface  charges  and  cur¬ 
rents  is  Love's  equivalence  principle  [24]  .  It  is  a  special  case  of  the 
general  field  equivalence  principle  [25]  in  that  a  null  field  is  postuated 

inside  the  surface  S,  and  the  charges  and  currents  of  Equations  (19 )-(22) 

21 


are  required  on  S  to  maintain  the  original  field  E,  H  outside  S.  Hone'- 
an  application  of  the  Huygens-Fresnel  integrals  using  these  surface  cu; 
rents  to  find  the  fields  at  a  point  inside  S  would  necessarily  yield  a 
null  result.  This  observation  is  important  when  applying  the  Huygens- 
Fresnel  integrals  to  find  specially  defined  fields  inside  the  rad one  as 
discussed  below. 

The  integral  representations  in  Equations  (14)  and  (15)  for  the 

fields  radiated  by  Antenna  "a"  into  the  homogeneous,  unbounded  medium 

outside  the  radome  are  exact  for  points  ui  to  and  including  the  surface 

S;  however,  the  difficulty  in  analysis  is  the  accurate  determination  of 

fields  E,  H  on  S .  To  gain  further  insight  into  this  problem,  consider 

the  geometry  of  Figure  4.  Define  primary  and  scattered  fields  of  the 

antenna  with  respect  to  surfaces  S  and  S2  as  was  done  in  Equations  (6) 

and  (7)  of  the  previous  section.  Let  the  primary  fields  E  ,  H  on  S, 

— ao  -ao  2 

be  known.  Then  the  primary  fields  incident  on  surface  S  can  be  found 
via  Equations  (14)  and  (15)  above  where  the  surface  of  integration  is 
now  S„  of  Figure  4;  i.e.. 


E 

—ao 


(S1> 


_1_ 
4  T 


v„. 

b2 


E  , 
— ao 


H  }dS 
— ao 


(26 


where  the  Huygens-Fresnel  operator  {  }  is  defined  here  to  be  the  in¬ 
tegrand  of  Equation  (14).  A  similar  expression  holds  for  the  magnetic 

field  intensity  H  .  By  definition,  the  fields  E  ,  H  are  those  that 
— ao  — ao  — ao 

would  exist  at  in  the  absence  of  the  dielectric. 

The  first-order  scattered  fields  E  ,  H  on  S  are  those  caused 

—as  —as  1 

by  the  boundary  and  the  media  outside  S  .  A  very  common  approxima¬ 
tion  for  this  complicated  scattering  process  utilizes  the  results 


of  plane  wave  interaction  with  a  dielectric  interface:  at  each  point  on 


S1,  the  incident  fields  E^q,  ate  approximated  locally  as  a  plane  wave 

and  the  dielectric  is  assumed  to  have  the  same  effect  as  a  plane  dielec¬ 
tric  panel.  The  normal  to  the  dielectric  surface  and  the  direction  of 
propagation  of  the  incident  fields  define  the  plane  of  incidence  and  the 
angle  of  incidence.  The  incident  electric  field  is  resolved  into  compo¬ 
nents  perpendicular  to  and  parallel  to  the  plane  of  incidence  as  illus¬ 
trated  in  Figure  5.  The  flat  panel  reflection  coefficients  R  ,  R  are 

1  il 

then  applied  to  the  incident  fields  to  find  estimates  of  the  first-order 
scattered  fields;  i.e. 


E 

—as 


a  (E 
1  — ao 


ai)Rx 


+  a  . 


(E 
— ao 


ii)  R  I.  =  E 

-ao 


(27) 


where  the  dyadic  R  is  defined. 


-  •  ;  na  *_  •. 


flat 

the 


panel  transmission  coefficients  T^,  T||  may  also  be  used 
fields  transmitted  to  the  outside  surface;  i.e., 


t 

E 

— ao 


a  (E 
l  — ao 


ai> 


T  +  a 
i 


,(E 

I  — ao 


a,|) 


A 

T„  =  E 
II  -ao 


(28) 


Ht 
— ao 


=  a  (H 
1  — ao 


V 


+  a 


i(H 

|  “-ao 


,,)T  =  H 

I  1  -ao 


(29) 


When  E^0>  H^q  are  used  in  Equations  (14)  and  (15) ,  a  first-order  approx¬ 
imation  to  the  radiated  fields  outside  the  dielectric  is  obtained. 

The  procedure  just  described  is  that  used  earlier  by  Paris  [13] 
and  emulated  by  Wu  and  Rudduck  [14]  and  Chesnut  (16) .  The  direction  of 
propagation  k  of  the  incident  field  was  taken  to  be  that  of  the  Poynting 


FIGURE  5.  PLANE  WAVE  PROPAGATION  THROUGH  AN 
INFINITE  PLANE  SHEET. 
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V 


vector  or  real  power  flow.  Also,  since  only  far-field  patterns  were  of 
interest,  the  usual  simplifying  approximations  were  applied  to  Equations 


(14)  and  (15)  to  obtain  the  following  far-field  approximations 


p  (RR  )  =  Z e"3kR 
hfiRRl]  4nR  6 


(Exn) xR 


[nxH  -  ( (nxH) *R^) R1  + 


dv  (30) 


Hff(RR1)  = 


R1  X  £ff 


(31) 


where  n  =  /p/e  is  the  characteristic  impedance  of  free  space  and  where 
o 

the  geometry  variables  are  defined  in  Figure  4.  Note  that  the  second  term 
in  the  integrand  of  Equation  (30)  is  the  radial  component  of  n  x  H  which 
is  subtracted  to  ensure  that  only  transverse  (to  R^)  components  of  the 
fields  are  found  in  the  far  zone.  Paris  used  the  fields  given  by  Equa¬ 
tions  (28)  and  (29)  in  the  integrand  of  Equation  (30) . 

The  first-order  scattered  fields  E  ,  H  on  surface  S,  also  con- 

— as  —as  1 

tribute  to  the  incident  fields  at  the  surface  as  indicated  by  the 
dashed  secondary  ray  in  Figure  4.  To  formulate  this  contribution,  con¬ 
sider  the  computation  of  the  field  E  at  an  interior  point  P  of  the 

3S 

volume  enclosed  by  S, .  Since  all  sources  of  E  ,  H  are  outside  S, , 

1  -as  -as  1 


E  (P)  =  j- 
— as  4  tt 


H  {-n,  E  ,  H  }dS 
F  —as  —as 


(32) 


where  -n  is  the  unit  inward  normal  to  and  T=l.  On  surface  write 


E  ,  H  as 
—as  —as 


E  (S)  =  E  •  R  +  E  (33) 

—as  1  — ao  =  —as 


H  (S,)  =  H  •  R'  +  H 
—as  1  — ao  =  —as 
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(34) 


where  R,  R'  are  flat  panel  reflection  dyadics  already  defined,  and  where 

II  II 

E  ,  H  represent  the  contributions  to  E  ,  H  at  each  point  on  S  from 
—as  —as  —as  —as  1 

radiation  by  the  fields  at  every  other  point.  Substitute  Equations  (33) 
and  (34)  into  Equation  (32)  with  T  =  2  for  P  on  S  to  yield 


E  (S. )  =  ~- 
— as  1  2it 


H_{-n,  E  ‘  R  +  E  ,  H  •  R'  +  H  }dS  (35) 

F  — ao  =  —as  — ao  =  —as 


Carry  out  the  integration  on  the  directly  reflected  terms  to  yield 


H_{-n,  E  ,  H  }dS  (36) 

,  F  —as  —as 

S1 

where  S|  is  surface  with  the  point  of  interest  excluded. 

II  M 

The  fields  E  ,  H  are,  of  course,  unknown;  however,  a  first-order 
—as  —as 

n 

approximation  to  E  at  a  point  P  on  S  is 
—as  1 

H  {-n ,  E  •  R,  H  •  R'}dS  (37) 

,  F  -ao  =  — ao  - 

1 

That  is,  the  directly  reflected  fields  at  every  other  point  on  are  used 

•I  II  II 

to  determine  E  at  the  point  of  interest  on  S, .  The  fields  E  ,  H  at 
-as  1  -as  -as 

each  point  on  S,  should  then  be  added  to  incident  fields  E  ,  H  to  pro- 
1  -ao  -ao 

duce  a  second  order  estimate  to  the  fields  on  the  outside  surface  of  the 
dielectric  via  Equations  (28)  and  (29).  Hence,  an  approximate,  iterative 
procedure  to  compute  the  effects  of  the  first-order  scattered  fields  on 
the  final  radiated  fields  is  formulated. 

An  important  hypothesis  for  the  validity  of  Equations  (14)  and 
(15)  is  that  the  fields  E,  H  on  S  satisfy  Maxwell's  equations,  including 
the  equations  of  continuity.  If  this  hypothesis  is  not  met,  the  fields 


E  (S  )  “ 
—as  1  2tt 


E  (S  )  =  E 
—as  1  — ao 


R  +  ~~ 

=  2  IT 


r - 


computed  at  points  not  on  S  will  still  satisfy  Maxwell's  equations  pro- 
f  vided  only  that  the  equations  of  continuity  are  satisfied  on  S;  e.g.,  if 

j  surface  currents  are  terminated  abruptly,  then  a  line  of  surface  charge 

i 

(  must  be  explicitly  included  in  the  field  expressions  [26] .  If  the  hy¬ 

pothesis  that  E,  H  on  S  satisfy  Maxwell’s  equations  is  met,  then  the  equa¬ 
tions  of  continuity  are  automatically  satisfied;  moreover,  the  expressions 
in  Equations  (14)  and  (15)  are  valid  for  points  on  surface  S  as  well  as  for 
points  outside.  Now,  the  fields  given  by  Equations  (26; ,  (33),  (34),  and 

(37)  do  satisfy  Maxwell's  equations;  however,  when  the  approximation  in 
Equation  (37)  is  substituted  into  Equation  (33),  the  latter  may  no  longer 
satisfy  Maxwell's  equations  on  s  because  of  the  first  term.  Hence,  a 

number  of  iterations  may  be  required  before  the  trie  values  of  E  ,  H 

— a  s  —a  ? 

on  are  found. 

Other  approaches  to  the  determination  of  the  first-order  scattered 

fields  E  ,  H  include  the  integral  equation  formulation  discussed  by 
sis  as 

Poggio  and  Miller  [22]  for  scattering  from  dielectric  bodies.  In  this 
approach  for  the  geometry  of  Figure  4,  integral  expressions  are  written 
in  each  region  for  the  electric  and  magnetic  fields  in  terms  of  the  equi¬ 
valent  surface  currents  and  charges  via  Equations  (14)  and  (15)  .  In 
Region  I  there  are  added  terms  TE^ ,  TH  in  the  integral  expressions  due 
to  the  incident  fields.  The  boundary  conditions  at  each  material  inter¬ 
face  are  then  written  using  the  integral  expressions,  resulting  in  a 
system  of  integral  equations  in  the  unknown  surface  charges  and  currents. 
In  principle,  the  equations  can  be  solved  by  the  method  of  moments  [27] 
or  by  iterative  methods  such  as  discussed  above.  The  radiation  fields 
in  Region  III  could  then  be  found  directly  from  the  solutions  for  the 
surface  currents  and  charges  on  via  Equations  (14)  and  (15) ;  however, 
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the  contributions  of  the  higher  order  scattered  fields  woulu  not  be 
included. 

The  techniques  discussed  in  this  section  for  the  determination  of 

the  fields  on  a  specified  surface  are  directly  applicable  to  the  receiving 

formulation  of  the  previous  section.  Most  notably,  in  Equation  (12),  it 

is  necessary  to  find  the  fields  e5  ,  h5  of  Antenna  ”b”  on  the  surface 

-bo  -bo 

S  which  coincides  with  the  radiating  aperture  of  Antenna  "a".  When 
ap 

Antenna  "b"  is  greatly  removed,  the  fields  incident  on  the  outside  of  the 

dielectric  are  those  of  an  infinite  plane  wave.  Appropriate  adaptations 

of  Equations  (28)  and  (29)  could  be  used  to  estimate  the  fields  E^,  H^Q 

on  the  inner  surface  of  the  dielectric  enclosure.  The  Huygens-Fresnel 

integrals  of  Equations  (14)  and  (15)  could  then  be  applied  to  transform 

these  fields  to  the  surface  S  .  This  approach  is  essentially  that  deve- 

ap 

loped  by  Tricoles  [28] .  An  integral  equation  approach  could  also  be  used, 
including  the  one  described  by  Van  Doeren  [10] . 

XV.  EQUIVALENCE  OF  FORMULATIONS 

It  is  well  known  that  the  receiving  and  transmitting  patterns  of 
an  antenna  in  free  space  are  identical.  In  this  section,  it  is  theoreti¬ 
cally  demonstrated  that  the  receiving  and  transmitting  patterns  for  an 
antenna  enclosed  by  a  dielectric  radome  as  in  Figure  1  are  also  identical, 
and  that  the  receiving  and  transmitting  formulations  of  radome  analysis 
lead  to  identical  results.  Establishment  of  the  equivalence  yields  in¬ 
sight  into  the  computational  advantages  of  one  formulation  over  the  other 
and  provides  a  means  to  compare  numerical  results  obtained  using  the  two 
approaches . 
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Consider  the  geometry  of  Figure  6  where  n  represents  an  inf. in' t<- 

b 

simal  electric  current  source  of  strength  I,  ,  lerioth  t f ,  and  or  ienta*  3 or, 

D 

located  at  the  r'oint  ( x',y’,z ')  which  is  a  great  distance  P.  fror  t 

origin  of  coordinates.  For  definiteness,  n,  is  transverse  to  R  ,  win  vh 
'  b  1 

is  the  orientation  that  would  normally  be  assigned  for  the  measurer'  :it 

of  the  radiation  pattern  of  Antenna  "a”  enclosed  by  the  surface  S0.  The 

fields  E,  H  at  (x^y'.z1)  due  to  Antenna  "a"  are  given  by  Equations  (14' 

and  (15) .  For  clarity,  let  E^,  represent  the  fields  oi:  produced 

by  Antenna  "a"  when  it  is  transmitting. 

Let  E  ,  H  be  the  fields  produced  on  S  by  the  current  element 
R  R  3 

n  .  An  application  of  the  general  reciprocity  theorem  to  the  unbounded 
b 

region  yields 


(E  x  ~  E  *  H  )  •  n  dS  = 

-T  — R  — R  -T 


n,  I  LI 
b  b 


A  second  application  of  the  reciprocity  theorem  to  Antenna  "a”  inside? 

the  surface  as  was  done  in  the  discussion  leading  up  to  Equation  ( 1 3) 

yields  the  result  that  the  complex  voltage  response  of  Antenna  "a"  to 

the  fields  E  ,  H  incident  on  S  is  proportional  to  the  integral  in 
R  K  3 

Equation  (38) ;  hence. 


V  -  X  Li.  E  •  n  =  C 
R  b  —  b 


where  C  is  a  complex  constant. 


(L  x  •  L  x  L ) 

-T  -R  — R  -T 


•  n  dS 


Equation  (39)  provides  the  connection  between  the  receiving  and 
transmitting  formulations  and,  as  such,  represents  the  third  facet  of 
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FIGURE  6.  GEOMETRY  USED  TO  ESTABLISH  EQUIVALENCE 
BETWEEN  THE  RECEIVING  AND  TRANSMITTING 
FORMULATIONS  OF  RADOME  ANALYSIS. 
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the  theory  of  radome  analysis.  Basically,  it  states  that  if  E,  H  are  tin' 
fields  of  Antenna  "a"  at  (x',y',z')  as  determined  via  the  Huygens-Fresne) 
principle,  then  the  voltage  response  of  Antenna  "a"  (when  receiving)  to 
an  elementary  current  source  at  (x’,y',z')  is  given  succinctly  by  Equa¬ 
tion  (39)  as  1  At  E  •  n,  .  An  antenna  more  complicated  than  the-  eleir.cn- 
b  —  b 

tary  current  source  could  be  assumed  for  Antenna  "b":  for  the  far-zone 
case,  the  same  results  would  be  obtained  since  the  radiation  and  recep¬ 
tion  properties  of  the  antenna  can  be  embodied  in  a  complex  effective 

vector  length  h  which  is  entirely  analogous  to  n  [29] ;  for  the  near- 

fa 

field  case,  the  approach  described  by  Paris  et  al.  [30]  would  be  required 
Even  though  Equation  (39)  is  valid  regardless  of  the  separation 
and  orientation  of  the  elementary  source,  it  is  desired  here  to  explicit¬ 
ly  establish  the  second  equality  in  that  equation  for  the  far-zone  case. 

To  this  end,  it  is  noted  that  the  fields  E1,  H1  on  S,  are  those  of  a 

HR  — R  3 

plane  electromagnetic  wave  propagating  in  the  k  =  -R^  direction  with  po¬ 
larization  properties  given  by  n,  according  to 

b 


E1 

-R 


jwu 


e jkk  •  o 


(4  ' 


i 

-R 


k  x  E 


o 


jkk  •  p 
e  -  - 


V 


(41) 


where  k  =  wu/n  =  2n/X  and  where  p  is  the  position  vector  from  the  origin 
o 

of  coordinates  to  the  source  point  on  as  indicated  in  Figure  C.  The 
constant  C  in  Equation  (39)  will  be  chosen  to  conform  to  these  explicit 
expressions  for  the  incident  plane  wave.  When  Equation  (14)  is  substi¬ 
tuted  for  E  in  I  A£  E  •  n  ,  there  results 
~  b  —  d 
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IAH 

j  = 

R  4ir 


-jkR  .  ~ 

.  .  e  ikk-p 

[-DWU  — - —  e 


£  %  •  <n  *  v  +  v  (<n  x  v 


-gkR 

}ke  ]kk • p , 

x  R,  — - -  e  -  +  n 


1  R 


)+  nb  •  Vi|/ff  (n*^)  JdS 


(42) 


where  the  following  asymptotic  relations  have  been  used: 


r  ~  R  -  R  *p 


(43) 


-]kr  -}kR 
e  e  ikk*P 

- e 


(44) 


„ ,  „ ,  '  e  ■'kR  jkk*p 

Ey  ~  — - —  e 


(45) 


Applying  the  vector  identity  a*bx£=b*£xa^to  the  first 
term  in  the  integrand  of  Equation  (42)  yields 


-jkR  '  „  -jkR  *  „  .  “ 

e  lkk’f.  ,  .  e  lkk-p  ,  , , 

-iwp  — - —  e  -  nb  •  (n  x  H^)  =  -jwy  — - —  n  -  x  e  -  (46) 


E  x  H  •  n 
-R  -T 


(47) 


i.e.,  it  is  equal  to  the  negative  of  the  second  term  of  the  integrand  in 
Equation  (39) .  Repeated  application  of  the  same  vector  identity  to  the 
second  term  in  the  integrand  of  Equation  (42)  yields 


n  •(<»  x  E_) x  R  x<jk  ^  R 

b  — T  1  R  '  -T  R  1 


x  V' 


(48) 


El  x  H 
-T  -R 


(49) 


1 3 


I 


i.e.,  it  is  equal  to  the  negative  of  the  first  term  in  the  integrand 

of  Equation  (39) .  The  third  term  in  the  integrand  of  Equation  (42)  is 

identically  zero  since  n,  is  perpendicular  to  Vip  .  Collecting  results 

b  —  It 

yields  the  final  desired  result  where  the  complex  constant  C  is  chosen 
to  be 


C 


I  Ax. 
b 

4n 


(50) 


to  ensure  equality. 


A  concise  application  of  the  foregoing  analysis  is  to  show  that 

the  response  of  Antenna  "a"  to  the  incident  plane  wave  depends  only  on 

the  incident  fields  E1,  H1  and  not  on  the  scattered  fields  E  ,  H  . 

R  R  RS  KS 

Direct  examination  of  the  first  equality  in  Equation  (39),  as  well  as 
examination  of  Equations  (47)  and  (49) ,  revelas  that  this  is  indeed  the 
case.  The  total  fields  E  ,  H  or.  S  are  given  by  the  superposition  of 

r\  i\  j 

the  incident  and  the  scattered  fields  as 


E  =  E  +  E 
— R  ~P.  -RS 


(51) 


H  =  H  +  H  _ 
-R  -R  -RS 


(52) 


Substituting  into  the  reciprocity  integral  of  Equation  (39)  yields 


VR  "  C 


(E„  x  -  E^  x  H) -n  dS  +  C  j  (E„  x  H  -  E  x  H  ) -n  as (53) 
-T  -R  — R  -T  -T  -RS  — RS  -T 

s  J  s 

3  3 


The  second  integral  must  be  identically  zero  because  of  the  equivalence 


established  explicitly  above;  also,  it  must  be  zero  by  virtue  of  the 


reciprocity  theorein  itself  as  discussed  in  Section  II  above  following 
Equation  (10) . 

V.  CONCLUSIONS 

The  theory  of  radome  analysis  is  based  entirely  on  the  recipro¬ 
city  theorem  and  the  Huygens-Fresnel  (or  Kirchof f-Huygens)  principle, 
both  of  which  are  derivable  from  Maxwell’s  equations  via  a  vector  Green's 
theorem  and  the  divergence  theorein  (Gauss'  Law).  All  methods  of  radome 
analysis  can  be  cast  in  terms  of  the  theory  presented  to  provide  a  common 
basis  of  understanding,  to  clarify  any  approximations,  and  to  provide  a 
basis  of  comparison  for  the  procedures  used  and  numerical  results  obtained. 

The  theory  presented  provides  the  correct  framework  in  which  to 
think  about  radome  analysis.  Its  understanding  is  of  paramount  impor¬ 
tance  to  the  correct  development  of  new  analysis  methods  which  may  util¬ 
ize  modern  numerical  methods  such  as  method  of  moments,  (7TD,  and  hybrids 
of  the  two.  The  theory  provides  the  basis  for  the  development  of  even 
newer  methods  of  analysis  which  may  rely  on  special  computer  architec¬ 
tures  or  networks.  Its  presentation  here  provides  a  solid  theoretical 
foundation  for  future  thrusts  in  this  important  area. 
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Abstract 

A  planar  aperture  synthesis  procedure  which  predicts  measured 
radiation  pattern  data  accurately  and  which  exploits  known  physical 
parameters  of  the  actual  antenna  is  described  for  use  in  computer- 
aided  radome  analysis.  The  plane  wave  spectrum  (PWS)  representation, 
with  the  geometrical  optics  approximation,  is  used  to  characterize 
radiation  from  the  array  of  four  identical,  discrete,  conical  horn 
elements  that  was  studied.  Solutions  for  the  PWS  of  each  element 
from  measured  array  patterns  over  the  visible  region  are  presented, 
and  a  digital  signal  processing  algorithm  is  described  for  extrapo¬ 
lating  the  aperture-limited  PWS  into  the  evanescent  region  as  re¬ 
quired  to  determine  the  near  field  of  each  element.  The  element 
near  fields,  having  bounded  support,  are  combined  to  produce  a  near 
field  for  the  complete  array.  The  array  near  field  is  used  in  a 
computer-aided  radome  analysis  to  demonstrate  the  accuracy  to  which 
the  measured  antenna  patterns  are  predicted  for  the  case  of  a  free 
space  radome.  For  completeness,  comparisons  of  measured  and  computed 
patterns  for  a  tangent  ogive  Rexolite  (e  =  2.54)  radome  are  presents 
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INTRODUCTION 


\ 


Valid  comparisons  of  the  accuracies  of  radome  analysis  methods 
using  measured  radiation  patterns  as  true  data  [1]  require  that  each 
method  under  consideration  be  capable  of  predicting  accurately  from  a 
suitable  antenna  near  field  the  antenna  patterns  measured  in  the  ab¬ 
sence  of  the  radome;  otherwise,  any  differences  between  measured  and 
calculated  patterns  with  the  radome  in  place  would  be  as  attributable 
to  the  inadequacy  of  the  antenna  model  as  to  the  radome  analysis 
method  itself. 

When  the  measured  amplitude  and  phase  of  the  antenna  fields 
are  known  over  a  complete  sphere  enclosing  the  antenna,  a  complete 
aperture  synthesis  method  such  as  that  described  by  Ludwig  can  be 
used  [2].  But  when,  as  in  many  cases  of  practical  importance,  the 
measured  pattern  data  is  limited  to  amplitude-only,  principal  plane 
patterns,  a  different  approach  and  one  which  exploits  all  available 
information  about  the  antenna  being  modelled  must  be  used. 

Such  a  method  of  planar  aperture  synthesis  for  the  four-horn 
monopulse  antenna  of  Figure  1  using  measured  principal  plane  ampli¬ 
tude  patterns  of  the  sum,  elevation  difference,  and  azimuth  differ¬ 
ence  channels  is  described  in  what  follows.  The  radius  (a  =  .74 A) 

and  element  half-spacing  (d  =  d  =  .95A)  of  each  vertically  (y) 

x  y 

polarized  conical  horn  element  (10°  flare  angle)  are  used  in  the 
antenna  model  to  describe  radiation  from  it  in  the  sinusoidal  steady 
state  (time  variations  of  the  form  e^ut  and  suppressed)  .  The  aper¬ 
ture  fields  of  each  identical  element  in  the  four-horn  array  are 
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1.  Geometry  of  Four— Horn  Monopulse  Antenna 


represented  by  their  plane  wave  spectra  [3]  as  determined  in  the  prin¬ 
cipal  planes  from  the  measured  patterns  and  the  above  physical  para¬ 
meters  of  the  actual  antenna.  The  two-dimensional  plane  wave  spectra 
in  the  other  radiating  portions  of  the  wavenumber  k^k^-plane  are  deter¬ 
mined  according  to  a  mathematical  model  deemed  characteristic  of  the 
actual  element  and  which  utilizes  principal  plane  spectra  only;  e.g., 
separable  spectra,  circularly  symmetric  spectra,  etc.  The  element  plane 
wave  spectra  in  the  non-radiating  or  evanescent  region  of  the  wavenumber 
plane  are  found  by  extrapolating  the  aperture-limited  spectra  from  the 
visible  to  the  invisible  region  using  an  algorithm  described  by  Papoulis 
for  one-dimensional  band-limited  functions  [4] .  The  element  plane  wave 
spectra  so  found  are  then  recombined  with  the  array  factor  to  produce 
ti.c  spectra  for  the  complete  antenna  which,  upon  Fourier  transformation, 
yields  a  near  field  for  the  antenna  which  accomplishes  the  desired  result. 

This  antenna  synthesis  problem  is  motivated  by  an  on-going  parame¬ 
tric  investigation  of  radome  analysis  methods  [1].  A  general  theory  has 
been  developed,  based  on  the  Huygens-Fresnel  principle  and  Lorentz  reci¬ 
procity,  which  embodies  the  various  methods  of  radome  analysis  [e.g.,  5-10]. 
Computer  codes  for  three  typical  methods  have  been  implemented.  Measured 
pattern  data  and  boresight  error  data  for  three  antenna  sizes  (small, 
medium,  and  large)  combined  in  fifteen  combinations  with  five  radomes  have 
been  obtained  to  serve  as  true  data  in  determining  the  accuracies  and 
ranges  of  validity  of  the  three  methods  of  analysis.  The  particular  mono¬ 
pulse  antenna  configuration  [11]  in  Figure  1  was  chosen  for  its  ruggedness, 
ease  of  fabrication,  and  ease  of  duplication  in  the  different  physical 
sizes.  Since  the  perturbations  in  the  patterns  and  boresight  errors  caused 
by  the  radome  are  small,  it  is  essential  to  the  success  of  the  research 


that  the  actual  antennas  be  accurately  modelled  in  the  analyses  so  that 


valid  comparisons  of  the  methods  can  be  made.  The  synthesis  procedure  was 
developed  specifically  for  this  purpose  and  for  this  antenna  configuration, 
but  it  has  wider  application  and  is  applicable  to  more  general  configurations 

The  plane  wave  spectrum  representation  of  the  antenna  fields  is  well 
suited  to  the  radome  analysis  application  because  of  the  computational  effi¬ 
ciency  which  can  accrue  as  a  result  of  the  Fourier  transform  relation!  . 
between  the  aperture  fields  and  the  plane  wave  spectra/far  fields.  And 
because  of  the  asympototic  relationships  between  the  far  fields  and  the 
plane  wave  spectra,  the  latter  can  be  mathematically  determined  from  the 
former,  at  least  over  the  visible  region  corresponding  to  the  hemisphere 
z>0.  But  before  the  desired  near  fields  can  be  determined,  it  is  neces¬ 
sary  to  assign  values  to  the  spectra  corresponding  to  the  evanescent  or 
non-radiating  modes.  Failure  to  do  so  may  result  in  a  near  field  which, 
because  it  contains  only  radiating  modes,  may  be  too  spread  out  to  fit 
within  the  confines  of  the  radome  for  analysis  purposes. 

In  the  synthesis  procedure  here,  values  are  assigned  to  the  eva¬ 
nescent  modes  by  extrapolating  the  plane  wave  spectra  from  the  known 
visible  region  into  the  evanescent  region.  Such  extrapolation  is  valid 
only  for  an  aperture-limited  spectrum;  i.e.,  a  spectrum  which  is  the 
Fourier  transform  (or  inverse  transform)  of  an  aperture  field  which  has 
bounded  support.  Such  a  spectrum  is  an  entire  function  [12]  of  the 


(complex)  wavenumbers  k^,k^  whose  value  in  one  region  of  the  complex 
k^  and  k^  planes  can,  in  principle,  be  determined  by  analytic  contin¬ 
uation  from  a  known  region  [13] .  The  algorithm  described  by  Papoulis 


provides  a  practical,  computer-aided  procedure  for  implementing  the 


extrapolation;  it  is  also  computationally  efficient  since  the  Fast 


Fourier  Transform  (FFT)  algorithm  can  be  exploited. 


This  synthesis  problem  is  actually  an  inverse  source  problem  such 
as  that  described  by  Schmidt-Weinmar  and  Baltes  [14]  in  which  prior  know¬ 
ledge  about  the  source  plays  an  important  role.  It  is  well  known  that 
such  problems  involve  the  mathematical  questions  of  existence,  uniqueness, 
and  stability  of  the  solution.  Small  errors  in  the  experimental  data, 
(errors  and  noise  are  inevitably  present)  can  lead  to  large  errors  in  the 
solution  unless  suitable  stabilizing  constraints  are  imposed;  i.e.,  un¬ 
less  additional  prior  knowledge  can  be  taken  for  granted  or  known  to  be 
fact.  The  procedure  described  here  utilizes  a  novel  combination  of  known 
and  assumed  data,  the  sum  of  which  provides  enough  information  to  produce 
a  stable  solution  to  the  inverse  problem  which  fits  all  the  prescribed 
known  data  —  but  which  may  not  be  unique.  The  question  of  uniqueness 
and  the  amount  of  information,  both  prior  knowledge  and  measured  data, 
needed  to  ensure  uniqueness  is  the  subject  of  a  current  investigation. 

This  paper  presents  the  solution  for  the  plane  wave  spectra  of 
the  antenna  in  terms  of  the  measured  patterns  and  the  physical  para¬ 
meters  of  the  antenna.  Spectrum  functions  separable  in  rectangular  and 
in  polar  coordinates  are  presented  which  permit  the  specification  of  the 
complete  radiating  spectrum  from  knowledge  of  it  in  only  the  principal 
planes.  The  extrapolation  procedure  in  both  one  dimension  and  in  two 
dimensions  is  detailed  and  used  to  effect  solutions  for  a  theoretical 
antenna,  whose  solution  is  known  exactly  for  verification  purposes, 
and  for  the  actual  antenna.  Computer  generated  results  are  presented 
in  the  forms  of  radiation  pattern  comparisons  and  three-dimensional 
plots  of  the  near  fields  and  plane  wave  spectra.  To  demonstrate  appli¬ 
cability  and  for  completeness,  some  radome  analysis  results  are  pre¬ 
sented  . 


45 


THEORETICAL  DERIVATIONS 


The  synthesis  procedure  is  carried  out  independently  in  t:.«- 

principal  p>lanes  to  determine  the  plane  wave  spectra  A  (k  ,  k  '  ,  A  , 

xe  x  y  ye 

and  aperture  fields  E  (x,y) ,  E  of  each  identical  element  in  the  P.  j> - 

xe  ye 

horn  array  of  Figure  1.  The  general  expression  for  the  radiation  f :  ■  :  1 
of  the  array  may  be  written  as  a  product  of  the  element  field  and  the 
array  factor  as 


Ehk  ,k  )  =  E  (k  ,  k  )  [a  e 
-ff  x  y  — eff  x  y  1 


jk  (d  k  +  d  k  ) 
o  x  x  y  y 


+  a2  b 


’k  (-d  k 
c  xx 


+  d  k 


jk  (-d  k  -  dk)  ]k  (d  k  -  d  k  i  1  (? 

,  o  xx  y  y  o  x  x  y  v 

+  a  e  +  a  e  1  • 

3  4 


where  k  =2r/l,  j  =  /-T,  the  a. ’s  are  the  complex  excitation  coefficients 
o  1 

of  the  elements  numbered  and  spaced  as  shown  in  Figure  1,  E  is  the 

— -e  f  f 

radiation  (electric)  field  of  the  element,  and  the  normalize^  wavenumber 


k  ,  k  are  related  to  spherical  angles  6,  $>  of  Figure  2  by 
x  y 


k  =  sinP  cosd 

x 


k  =  sinO  sinf  ,  f  3 

y 

at  least  for  the  visible  region  defined  bv  y'k  ‘+k  ‘  <  1 

•  x  y 

The  three  monopulse  channels  of  the  antenna  are  formed  by  phasir.c; 
the  elements  appropriately.  In  the  ideal  case,  |a  [=1  for  i=l-4  .  For 
the  sum  (I)  channel,  all  elements  are  assigned  equal  phase.  For  the  ele 
vation  difference  (A^)  channel,  a  =a2=-a  =-a  .  For  the  azimuth  differ¬ 
ence  ( )  channel,  a^=a^=-a^=-a^ .  For  perfect  assignment  of  the  exci- 


FIGURE  2.  COORDINATE  SYSTEM  FOR  ANTENNA  PATTERN 
MEASUREMENTS  ON  VERTICALLY  POLARIZED 
ANTENNA. 
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E  -  E  (k  ,k  ) (4  cot  k  d  k  cos  k  d  k 
—ft.  -t.fr  x  y  o  x  x  o 


=E,,(k,k)fj4ccsk  d  k  sir.  k 

T  ’  — F  •  f  I  X  V  -v  -v 


F  . , .  =  E  ,,(k  ,k  )  (-  4  sir.  k  d  k  cos  k  u  k  ) 

Ac  -off  x  y  o  x  x  c  v 


'■  r.or.- ideal  case  encountered  in  practice,  the  t-xci  tat  :  c  r.  ct-*f  f 

-it:.:,  values  may  vary  somewhat  from  the  ideal  value;  sc  that  th- 
■  !■  x  array  factor  of  Equation  (1)  must  be  used  in  sol  vine  for  tie  t 
:  field  and  plane  wave  spectra  of  measured  sum  and  different--  :  ~ 

The  radiation  field  of  the-  element  is  related  to  its  y  lari'-  ■■> 
spectra  according  to  [15] 


E  „(k  k  )  = 


x  [  ( 1+k  -k  )  A  -  k  k  A  ] 
z  x  xe  x  y  ye 


+  y[-k  k  a  +  (1+k  -k  )  a  ! 
x  y  xe  z  y  ye 


+  z[-k  (1+k  )A  -  k  (1+k  )  A  1 
x  z  xe  y  z  ye 


where  k  =Vl-k  -k  .  The  geometrical  optics  approximation  [If]  ha: 
z  x  ^ 

assumed  for  the  aperture  fields  .in  the  conical  horn  element;  viz., 


— ap  n 


*■.!  v‘  x,  y,  z  are  unit  vectors  in  the  rectangular  coordinate  rp  * . 
i'l  r-'  ,  ar.d  n  is  the  wave  impedance.  In  Equation  (T)  ,  the  rad  i  a 


of  the  forr  p  *  /r  has  been  suppressed  where  t  1  s  t  he  r.i  ii 


V 


distance  from  the  origin  of  antenna  coordinates  to  the  far-field  mea¬ 
surement  sphere  of  Figure  2. 

Equation  (7)  may  also  be  written  in  terms  of  the  transverse- 
spherical  comp’oner.ts  as 

( 1  +k  * 


-e  f  f 


(k  ,k  ) 
x  y 


- - - 1  v  ( k  A  +  k  A  )  +  $(-k  A  +  k  A  )  ] 

I — - —  x  xe  v  ye  y  xe  x  ye 

V'k  Z  +K 


('. J) 


x  y 

which  is  a  form  more  suitable 
sured  transverse  components 
V_  defined  by 


to  solution  for  A  ,A  in  terms  of  the  mc-a- 
xe  ye 

,  E  or  the  related  measured  quantities  V. , 


V.  = 
o 


k  >  0 

x 

k  <  0 

x 


V 

<? 


k 

x 

k 

x 


0 

0 


(11) 


For  k  =0  ( E-viane ) ,  the  definitions  above  apply  when  k  is  substituted 
x  -  y 

for  k  .  Let  AF  ,  i=l, 2, 3,  denote  the  array  factors  associated  with  the 
x  l 


..  ,  L  ,  L  channels,  respectively,  modelled  in  Equation 
solutions  for  the  element  spectra  are  given  by 


Then  the 


A  (k  ,k  ) 
xe  x  y 


(1+k  ) A  2+k  2  AF . 
z  x  y  i 


(12) 


=  -X 


k  V,  +  k  V.  . 
°i  x  $1 


(1+k  )4  "2~+k  zhF 

z  x  y  i 


A  (k  ,k  ) 

ye  x  y 


(13) 


for  all  k  ,k  sue!',  that  Oik  <1.  At  the  pole,  k  =1  and  the  spectra  mav 
x  y  z  z 

found  from  the  sum.  channel  measured  voltaaes  as 


A  (0,0) 
xe 


(0,0) 


"V(0'0) 


A  (0,0) 
ye 


V^1 (0,0) 

<pl 


vcx  (0,0) 


where  the  superscripts  denote  the  conventional  E  and  H  planes  [17] .  T! 
general  expressions  in  Equations  (12)  and  (13)  can  be  specialized  to  t! 
principal  planes  as  follows: 


H-Plane  (k  =0) : 

y 


A  (k  ,0)  = 
xe  x 


Ve-  (k  > 

01  X 


(l+/l-k  i)AF. 

X  1 


A  (k  ,0) 

ye  x 


v"(k  ) 

_ 4>i  x _ 

(l+/l-k  2)AF. 

X  1 


E-Plane  (k  =0)  : 
x 


A  (0,k  )  = 
xe  y 


“V  (k  } 

(1+A-k  2)  AF. 
y  i 


A  (0,k  )  = 
ye  y 


Vp. (k  ) 

ci  y 


(1+vl-ky^) AF 


When  numerically  evaluating  Equations  (1C) -(19) ,  the  channel 
dex  i  is  chosen  on  the  basis  of  h_ghest  signal  level.  For  example, 

1  channel  (i=l)  provides  the  best  a  .ta  for  small  k  in  the  H-;  lar,<. 


tn 


nowcver,  as 


k  i  increases,  the  sum,  channel  ami  litude  decreases  towa 
x 


first  null  in  the  array  factor  AF  while  the  channel  am:  1: 

-  AZ 


increases  toward  its  peak  at  k  ~  :sin  ( (4d  /a)  }  .  When  alter;. at: 

x  x 


between  the  data  sets,  it  is  important  to  ensure  that  the  measure  c.  d 
ference  channel  data  is  correctly  normalized  with  re-.:  >•.■'.  to  the  sum 
data  as  dictated  by  the  model  in  Equation  (1).  For  .  .  x^it a  tier, 

efficients,  the  ratio  of  difference  channel  resr-o:..--  ' 

sponse  is  given  in  each  principal  plane  by 
H-Plane : 

A_(k  ,0) 

-p— -  =  3  tan  k  d  k 
l.  (k  ,0)  o  x  x 


.  .  .a.L.  .  l:*  I  t* 


E-Planc  : 


A  (0,k  r) 

=  j  tan  k  d  k 
I  (C,k  )  J  o  y  v 


These  relationships  can  be  used  to  adjust  tin-  iuta  for  consis 

tent  solutions  of  the  element  spectra. 

The  measured  quantities  V,  ,  V  ■  . ons  \’6)-(19)  are  corri 

6i  :  . 

plex  quantities  usually  measured  us;:.  .  amplitude  receiver.  In 

those  myriad  cases  where  far-fieid  j  ;  ■  ;.ita  are  not  valid  due  to 

range  imperf ections ,  unstable  microwave  sources,  etc.,  some  assume -ion 
must  be  made  in  assigning  phase  vaiues  to  the  measured  amplitude  data 
each  angular  direction.  In  this  investigation,  phase  data  were  assign 
as  dictated  by  the  array  factor  for  each  channel.  The  phase  of  the 


(25)  depend  on  the  excitation  coefficients;  hence,  an  iterative  procedure 
is  indicated.  The  individual  excitation  coefficients  can  then  be  found, 
to  within  a  constant,  from  the  four  above  equations.  For  definiteness, 
a  value  of  (1+jO)  may  be  assigned  to  a  .  The  solutions  for  the  other 
coefficients  follow  as: 


given  by 


c  =  p  (0)  q  (0) 
l  x  x 
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C  =  p  (0)  q  (0) 
2  ry  y 


The  Fourier  transformation  above  implies  that  the  spectra  are  known  for 

all  k  (k  )  on  the  real  line;  however,  Equations  (16) -(19)  produce  solu- 
x  y 

tions  only  in  the  visible  region  |k^|  <  1,  ]k  |  <  1  of  the  principal  plants. 
The  extrapolation  technique  for  assigning  values  to  the  spectra  outside 
the  visible  region  will  be  deferred  until  after  the  discussion  immediately 
below  concerning  another  type  of  separability  for  spectra  and  near  fields 
which  will  be  applied  during  the  synthesis  procedure. 

The  TE  circular  waveguide  mode  (y-polarized)  provides  a  model 
for  an  element  near  field  that  is  separable  in  cylindrical  coordinates  c 
and  and  which  depends  only  on  knowledge  of  the  principal  plane  PWS . 

Assume  that  the  element  near  field  can  be  written  as 


E  (p,<p)  =  p  f  (p)  g  (4>)  +  <p  f  (p)  g  (<p) 
— nf  p  p  <t>  q 


(37) 


Convert  to  rectangular  coordinates  and  assume,  as  in  the  case  of  the  TF^. 
waveguide  mode,  that 


g  (4>)  =  sirnf  (?q' 

P 


g  ,  (<f>)  =  cos4>  (  3°) 
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The  element  (tangential)  near  field  can  then  be  expressed  in  rectangular 
coordinates  as 


„  f  -f  „  f  +fA  f  -f . 

E  (x, y)  =  x  P  ^  sin2$  +  y  (— — - - ^ -  cos2f) 

—nr  z  z  z 


where  x=pcos<(>,  y=psin<{)  as  usual.  Inverse  Fourier  transform  these  rec¬ 


tangular  components  to  yield  integrals 


A  (k  ,k  )  = 
x  x  y 


"p  $  _  ,  f  .  „ .  j 2irp  k  cos($-£) 
—  -x  Pdp  j  sm2i}>  e  p  < 


„  „  ,  .  p  ip  ,  j27tpk  cos  ^*5 L, 

A  (k  ,k  )  =  — — —  pdp  eJ  p  T  *'d$ 

yxyJ2  J 


fp  f(f)  ,  f  j2irp  k  cos(<J>-£)  ,  . 

— Pdp  cos2<£  e  p  d<J> 


Carry  out  the  integrations  in  4>  [19]  to  obtain 


A  (k  ,k  )  =  — 2tt  *  \  H  (k  ) 
x  x  y  k  0  P 


k  2-k  2 


A  (k  ,k  )  =  it  H  (k  )  +  tt 
y  x  y  Op 


— H  (k  ) 

2  2  p 


where  k  =  /k  z+k  *  and  where  the  Hankel  transforms  H  ,  are  defined 
p  x  y  0  2 

by  [20] 
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VV  i  <W  Jo{kpp)  pdp 


H  (k  )  A  (f  -f.)  pdp 

2  p  =  I  P  9  2  p 


Specialize  Equations  (43)  and  (44)  to  the  principal  planes,  and  solve  the 
resulting  system  of  equations  for  Hq,  to  yield 


VV  = 


A  (k  ,0)  +  A  (0,k  ) 

v  p  y  p 


H2(kP)  = 


A  (k  ,0)  -  A  (0  ,k  ) 

y  p _ y  P 


Thus,  the  two-dimensional  plane  wave  spectra  A  ,A  are  expressible  in 

x  y 

terns  of  only  the  principal  plane  spectra  via  Equations  (43) -(48). 

Since  the  assumed  trigonometric  variations  of  g  ,g  in  Equat’on= 

P  9 

(38)  and  (39)  do  not  produce  any  x-component  of  near  field  aiona  x=fi  or 
y=0,  a  similar  analysis  can  be  carried  out  for  an  orthogonal  TE-type  mode 
which  will  account  for  the  presence  of  such  cross-polarized  component?. 
Write  the  total  element  fields  as  the  sum  of  two  components 


E  .  =  p  (E  +E’)  +  <f>  (F,  +E;) 
— en  f  p  p  9  b 


where  E  ,  E  are  given  in  Equation  (37)  and  where 

P  9 


E'  =  f'(o)  g’ (9)  =  f'  cos9 
pop  p 


E'  =  f'(p)  a' It)  =  f'  sini 
9  9  9  9 


I 


Follow  through  with  the  analysis  and  define  Hankel  transforms  by 


G  (k  )  = 
0  P 


(f'-fM  J  (k  p)  pdp 
p  (J>  OP 


(52) 


G  (k  )  = 
2  p 


(f*+f*)  J. (k  p)  pdp 
P  <j>  2  p 


(53) 


Combine  all  results  for  the  two  orthogonal  modes  to  yield  the  following 

expressions  for  the  element  plane  wave  spectra  in  terms  of  the  principal 

H  £  H  £ 

plane  spectra  A  (k  ) ,  A  (k  ) ,  A  (k  ) ,  A  (k  ) : 

xe  p  xe  p  ye  p  ye  p 


A  (k  ,k  ) 
xe  x  y 


A  +A  k  z-k 
xe  xe  .  x  y 
— — — — —  +  - - 


-A  \  k  k 

xe  xe  +  (aE  _rH 

1  k  z  ye  ye 


(54) 


AH  +ae 


(55) 


The  above  two  equations  constitute  a  reasonable  model  for  the  plane 
wave  spectra  of  a  circular  element  such  as  that  shown  in  Figure  1.  The 
principal  plane  spectra  ir  these  equations  are  found  from  measured  prin¬ 
cipal  plane  patterns  via  Equations  (16)-(19).  Note  that  the  above  equa¬ 
tions  reduce  identically  to  the  measured  spectra  in  the  principal  planes. 
Note  also  that  the  model  is  based  on  the  characteristic  modes  of  a  cir¬ 
cular  aperture  geometry  and,  in  the  absence  of  complete  measured  data 
which  could  perhaps  determine  exactly  the  plane  wave  spectra  off  the 
principal  planes,  this  model  is  more  suitable  for  the  antenna  described 
than  is  the  model  of  Equations  ( 31 )  —  ( 3 2 )  which  assumes  rectangular 
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Strictly  speaking,  the  latter  would  more  suitable  {■  r  j.; 


arability. 

at:  a  rectangular  element  geometry,  but  does  find  app  1  i  cat  i  <•,  <  r.  *  i: 
xtrapolation  procedure  described  next . 


EXTRAPOLATION  PROCEDURE 

The  extrapolation  procedure  is  best  explained  in  one  dinensi on 
before  extending  it  to  the  two-dimensional  case.  The  objective'  of 
procedure  is  to  assign  values  to  the  evanescent  modes  in  the  plane'  wav 
spectrum  of  the  element  such  that  the  known  spectrum  in  the  visible  re¬ 
gion  is  preserved  and  such  that  the  near  field  of  the  element  is  re¬ 
stricted  to  a  finite  portion  of  the  xy-plane.  Results  are  presented  in 
this  section  which  demonstrate  the  algorithm  for  beth  an  ideal  case  and 
a  case  which  uses  measured  data. 

The  flow  chart  in  Figure  3  summarizes  the  extrapolation  procedur 
for  the  case  of  plane  wave  spectrum  A^Ck^)  and  its  corresponding  near 
field  E  (x)=F{A  (k  )}.  The  zeroeth  estimate  to  the  actual  spectrum  is 

X  XX 

the  portion  known  only  over  radiating  wavenumbers  as  defined  by 


A  (k  )  = 
XO  X 


A  (k  )  ,  Ik  !<  1 

|  X  X  X 


,  elsewhere 


The  known  length  2a  of  the  interval  over  which  the  near  field  is  non¬ 
zero,  and  a  lower  bound  1-e  on  the  fraction  of  enerqy  in  the  near  field 

\ 

which  must  be  contained  in  this  interval  for  an  acceptably solution , 


constitute  the  remaining  elements  of  input  data.  y 

The  first  estimate  E'  of  the  near  field  is  found  by  Fourier 

xo 

\ 

transformation  of  A  .  Due  to  the  lack  of  radiating  modes  in  the 
xo 

% 

spectra,  the  near  field  will  be  sjVead  out  in  x.  That  portion  lyinc 


I  o 

(INITIAL  ESTIMATE) 


|  AXH(kx)|  |  -  EXNF'lx, 


ESTIMATE 


EXNF(x).  =  EXNF'(x)  ■  JT 


AXH"  =  ,>-1  {  EXNF,  } 


AXH,  =  AXH0  +  AXH",  "LT 


FIGURE  3.  FLOW  CHART  FOR  EXTRAPOLATION  ALGORITHM. 
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outside  I  x  j  <a  is  then  truncated  so  that  E  =Ef  r.  (x)  results,  wh'-re 

xo  xo  a 

l-  (x)  is  the  rectangular  function  whose  value  is  unity  for  jxU'a  and 
a 

zero  elsewhere.  The  resulting  near  field  E  is  inverse  Fourier  t ran  s - 

xo 

formed  to  produce  the  spectrum  A*  (k  )  which  does  contain  non-zero  eva- 

xl  x 

Descent  modes  due  to  the  bounded  support  of  E  .  The  evanescent  modes 

xo 

of  A’,  are  added  to  the  radiating  modes  A  to  produce  the  first-order 
xl  xo 

s'.v ctrum  A  ,  =A  +A'  ,  p,  (k  )  where  p,  (k  )  is  the  complement  cf  the 
xl  xo  xl  1  X  lx 

rectangular  function  whose  value  is  unity  for  |xj>a  and  zero  for 
The  process  is  repeated  until  the  fraction  of  energy  in  the  near  field 
F  (x)  on  the  interval  jxjsa  exceeds  the  lower  bound  specified.  The 
EFT  algorithm  is  used  to  effect  the  Fourier  transformations. 

The  algorithm  was  first  tested  on  an  ideal  four-horn  monopulse 
array  (Figure  1)  where  a  TE  ^  circular  waveguide  mode  was  assumed  pre¬ 
sent  in  each  element.  The  mode  was  y-polarized  but  then  rotated  by 
a=5°  counterclockwise  so  that  a  small  x-component  of  radiated  field 
would  be  present.  The  far  fields  of  the  array  were  computed  at  2°  in¬ 
crements  in  the  principal  planes  for  all  three  channels  over  a  dynamic 

range  of  40  decibels.  These  data  were  used  as  the  "measured"  data  V  , 

H  H 

V, .  in  Equations  (16) -(19)  to  produce  element  spectra  A  (k  ) ,  A  (k  ) , 

®i  xe  x  ye  x 

E  E 

A  (k  ) ,  A  (k  )  over  the  visible  region  in  each  plane  at  2°  increments 
xe  y  ye  y 

in  polar  angle  6.  Since  the  discrete  Fourier  transform  relationship 
between  spectrum  and  near  field  requires  samples  of  the  spectra  at  equal 
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where  6=sin  k  ,  N  =  number  of  angular  samples  (91)  ,  and  6  =  angle 

x  max 

corresponding  to  the  Nth  sample  (90°) .  The  interpolated  spectra  were 

then  operated  on  in  turn  by  the  extrapolation  algorithm. 

Some  computed  results  are  shown  in  Figure  4  where  comparisons  of 

the  true  and  synthesized  near  fields  of  the  element  in  the  x=G  plane  arc- 

presented.  The  true  x  and  y  field  components  are  simply  those  of  the 

rotated  TE  mode.  After  twenty-one  iterations  of  the  extrapolation 

algorithm,  the  concentration  factors  for  the  x  and  y  components  were 
E  E  -4 

=  .013  and  =  .330  x  10  ,  respectively.  (In  the  y=0  plane, 

H  H 

slightly  better  results  were  obtained;  i.e.,  =  .0153,  e  =  .982  x 

-5 

10  ) .  The  half-length  of  the  interval  occupied  by  the  element  was 

a=.74136A,  The  algorithm  was  quite  successful  in  concentrating  the  near 
field  into  the  desired  interval  and  in  predicting  the  spatial  variation 
of  the  near  field  reasonably  well.  The  sample  spacing  in  the  near  field 
was  Ay=A/8. 

The  element  near  fields  in  Figure  4  were  (inverse)  Fourier  trans¬ 
formed  to  produce  the  corresponding  spectra  at  equal  wavenumber  over  the 

range  Ik  1 <k  =(2  Ay/A)  ^=4 .  The  spectra  were  then  Fourier  interpo- 
1  y 1  ymax 

lated,  to  produce  their  values  at  2°  increments  over  the  range  ]6;£90° 


The  interpolated  values  were  used  in  the  antenna  model  of  Equation  (1) 

to  produce  the  sum  channel  E-plane  patterns  shown  in  Figure  5.  Differ- 

jr 

ences  between  the  true  and  synthesized  primary  E  components  are  hardly 

0 
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FIGURE  4.  COMPARISONS  OF  TRUE  AND  SYNTHESIZED  ('!  NEAR  FIELDS  OF  THEORETICAL 
FOUR  HORN  MONOPULSE  ANTENNA  ELEMENT  (n  -  5°f. 


LEGEND 


FIGURE  5.  COMPARISONS  OF  TRUE  AND  SYNTHESIZED  (')  E-PLANE  SUM 

PATTERNS  OF  THEORETICAL  FOUR-HORN  MONOPULSE  ANTENNA 


discernable  except  at  very  wide  angles.  Good  agreement  is  obtain*-  :  f 
E 

the  cross  E  component  where  all  amnlitudc-  levels  are  less  than  o:- 

t *\ 
f 

Similar  results  were  obtained  in  the  H-plane  and  for  the  other  tv 

riels.  In  conclusion,  the  synthesis  procedure  works  quite  well  in  or.- 

dimension  for  theoretical  data  having  a  dynamic  range  of  4c  decih*-! r . 

The  one-dimensional  synthesis  procedure  was  next  tested  usir.e 

principal  plane  patterns  obtained  at  35  GHz  for  the  actual  antenna  of 

Figure  1.  Figure  6  shows  the  first  and  final  estimate?  of  the-  clf-r.o  t 

near  field  in  the  x=0  plane .  After  51  iterations  of  the  extrapolation 

E  E 

alaorithm,  the  E-plane  concentration  factors  were  c  =.069  and  <;'=.3V'X 

x  y 

-3  H  H 

10  .  (The  H-plane  concentration  factors  were  e  =.50F  and  e  =.494 

x  y 

_ 

10  .)  The  sample  spacing  was  Ay=.025>.,  yielding  a  maximum  wavenumber 

k  =20  for  N=256  samples  over  the  total  interval  of  6.4’.  Rc-ferrincr 
ymax 

to  Figure  6,  it  is  seen  that  the  algorithm  successfully  concentrated 
the  near  field  into  the  interval  [ yj <a=. 741361  as  desired. 

Figure  7  shows  comparisons  of  the  measured  patterns  with,  those 
synthesized  from  the  extrapolated  element  spectra  as  explained  above 
for  the  ideal  array.  The  agreement  is  quite  good  in  the  E-plane, 
Similar  results  were  obtained  for  the  H-plane  and  for  the  other  two 
channels.  In  conclusion,  the  synthesis  procedure  works  well  for  mea¬ 
sured,  one-dimensional  data  having  a  40-dB  dynamic  range. 

TWO  DIMENSIONAL  EXTRAPOLATION 

The  one  dimensional  results  could  be  used  through  Equations  ’ 
and  (32)  to  generate  two-dimensional  spectra  and  near  fields;  however . 
if  the  pfi-senarability  model  is  used  as  preferred,  then  the  extra polar 
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FIGURE  6  FIRST  AND  FINAL  ESTIMATES  OF  APERTURF  FIE  L D  OF 

ACTUAL  ELEMENT  IN  E  PLAN! 
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procedure  must  be  carried  out  in  two  dimensions.*  This  section  presents 

results  of  two-dimensional  synthesis  of  the  element  spectra  and  fields 

of  the  antenna  in  Figure  1  using  measured  principal  plane  patterns. 

Figure  8  shows  the  element  near  field  E  (x,y)  and  the  correspond- 

ye 

ing  plane  wave  spectrum  A  (k  ,k  )  obtained  using  the  2-D  extrapolation 

ye  x  y 

algorithm.  The  element  principal  plane  spectra  were  first  determined  via 

Equations  (16) -(19)  from  the  measured  data.  These  results  were  then  used 

in  Equations  (54)  and  (55)  to  generate  the  2-D  plane  wave  spectra  over 

2  2 

the  visible  region  k  +k^  il .  The  extrapolation  algorithm  was  then 

applied  to  each  2D  spectrum  in  turn  using  a  2-D  Fast  Fourier  Transform. 

After  eleven  iterations,  the  concentration  factors  were  £  =.065  and 

x 

£  =.028;  i.e,,  for  the  y-component  E  (x,y),  98.2%  of  the  near-field 
y  ye 

energy  was  concentrated  in  the  circular  element  region  of  radius  a=.74316A. 

A  square  array  of  64  x  64  samples  spaced  at  Ax=Ay= . 07031 A  and  resulting 

in  k  =k  =7.11  was  used  for  each  component.  An  exponential  window 
xmax  ymax 

was  used  to  limit  the  near  fields  to  the  circular  region. 

The  resulting  element  spectra  were  next  combined  via  Equations  (1) 
and  (9)  to  produce  the  2D  plane  wave  spectra  for  the  three  monopulse  chan¬ 
nels.  Each  spectra  was  then  Fourier  transformed  to  obtain  the  correspond¬ 
ing  near  fields.  Some  results  for  the  sum  channel  are  presented  in 

£ 

Figure  9.  The  FWS  A  (k  ,k  )  shown  in  Figure  9(a)  has  maximum  wavenumbers 

y  x  y 

k  =k  =3.56.  Although  no  extrapolation  was  applied  to  the  PWS  con- 
xmax  ymax 

taining  the  array  factor,  the  actual  measured  spectra  in  the  principal 
planes  were  inserted  directly  into  the  final  PWS  to  guarantee  close  agree¬ 
ment  with  measured  patterns  as  extracted  from  the  2-D  PWS. 

*One-dimensional  extrapolation  could  possibly  be  used  by  utilizing  the 
Hankel  transform  in  place  of  the  Fourier  transform. 
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FIGURE  8 


V 


Figure  9(b)  shows  the  y-component  of  the  synthesized  near  field 
of  the  four-horn  monopulse  antenna  at  64  x  64  samples  spaced  at  Ax=/,y- 
. 14045A.  The  near  fields  of  the  four  circular  elements  are  clearly 
visible,  and  the  vast  majority  of  the  near-field  energy  is  concentrated 
in  the  central  portion  of  the  array. 

When  selected  cuts  were  extracted  from  the  2-D  phase  wave  spectra 
of  the  sum  channel  and  used  in  Equation  (9) ,  the  patterns  of  Figure  1' 
resulted.  Agreement  between  the  measured  and  synthesized  E-plane  sum 
patterns  shown  in  Figure  10(a)  is  excellent.  But  perhaps  more  impor¬ 
tantly,  the  agreement  in  the  4>=4 5 °  plane  is  also  excellent  as  shown  in 
Figure  10(b)  —  facts  which  attest  to  the  apparent  correctaness  of  the 
underlying  assumptions  used  throughout  the  synthesis  procedure.  Similar 
.esults  were  obtained  for  the  other  planes  and  other  channels. 

RADOME  ANALYSIS  EXAMPLE 

The  motivation  for  the  antenna  synthesis  procedure  described  above 
is  to  ensure  the  valid  assessment  of  the  accuracies  of  various  computer- 
aided  radome  analysis  techniques  [1],  one  of  which  is  described  in  tio  •- 
section . 

A  computationally  fast  receiving  formulation  for  radome  analysis 

[9]  is  illustrated  in  Figure  11.  The  antenna  near  fields  E  . ,  H  ,  are 

— Ti  — Ti 

represented  by  a  uniform  grid  of  sample  points  on  a  finite  planar  aper¬ 
ture  surface  S  placed  on  or  just  in  front  of  the  radiatinc  portion  of 
ap 

the  actual  antenna.  The  near  fields  are  assumed  to  be  the  sarn^  as  these 
produced  by  the  antenna  when  transmitting  in  the  absence  of  the  radome. 

Rays  representing  the  incident  plane  wave  (target  return)  are  traced  from 
each  point  in  the  aperture  in  the  direction  k  to  the  radome  wall.  The 
electromagnetic  field  associated  with  each  ray  is  weighted  by  the  insert  ion 
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(a)  E -Plana  Pattarni. 


legend 


(b)  Pattarm  in  «  •  46°  Plant 


FIGURE  10  COMPARISONS  OF  MEASURED  AND  SYNTHESIZED  H  SUM  PATTERNS 
OF  ACTUAL  ANTENNA  USING  2-D  EXTRAPOLATION 
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transmission  coefficients  for  field  components  parallel  and  perpendicular 


to  the  plane  of  incidence  at  each  point  on  the  radorr.e  wail.  The  response 

V  of  the  antenna  is  obtained  by  summing  uj  the  contributions  of  the  re- 
R 

ceived  fields  E  ,  H_  as  specified  by  Loren tz  reciprocity;  viz., 

K  K 

-  c 

V„.  (k)  =  C  j  (E  x  H  -  E  x  H  . )  •  n  da  ( >  •  i 

Ri  J  — Ti  — R  — R  — Tl 

S 

at; 

The  subscript  i=l,2,3  specifies  the  antenna  near  fields  for  the  I'.,  A  ,  ana 

EL 

A  channels,  respectively.  Antenna  patterns  can  be  computed  by  control- 
ling  the  direction  of  arrival  and  polarization  of  the  incident  plane  wave . 
Boresight  errors  in  the  two  monopulse  planes  can  also  be  computed  for 
specified  antenna/radome  orientations  by  determining  the  direction  of  arri¬ 
val  which  produces  the  difference  pattern  nulls. 

The  64  x  64  sample  array  of  Figure  9(b)  representing  the  near  field 
of  the  antenna  was  actually  too  large  to  fit  into  the  10. 38A-diameter , 

10 . 38  A  length,  tangent  ogive,  Rexolite  (e^=2.54,  tan  6=.002)  radome  used 
in  the  analysis  and  measurement.  Also,  since  six  such  complex  arrays  are 
used  in  the  program  for  the  three  antenna  channels,  core  memory  storage  was 
a  consideration.  Consequently,  only  the  central  49  x  49  sample  points  of 
the  near  fields  were  actually  used  in  the  radome  computations. 

Comparisons  of  the  E-plane  A  patterns,  measured  without  the  ra- 

ELi 

dome  and  computed  using  a  free  space  radome,  are  shown  in  Figure  12  to 
demonstrate  the  adequacy  of  the  synthesis  technique  and  the  small  effects 
of  the  near-field  truncation.  The  pattern  in  Figure  12  was  computed  as 
the  response  of  the  antenna  to  plane  waves  arriving  from  65  directions  in 
the  E-plane,  equally  spaced  in  k  =sin0  over  the  visible  region  |k  |sl.  The 
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FIGURE  12.  COMPARISONS  OF  MEASURED  AND  COMPUTED  (')  E-PLANE  DIFFERENCE 
PATTERNS  FOR  FREE  SPACE  RADOME. 
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PATTERNS  FOR  REXOLITE  P.ADOME. 


computed  patterns  were  then  Fourier  interpolated  via  Equation  (IB)  to  : r 
duce  sample  points  at  2°  increments  over  the  angular  range  '/  ’  *"9''°. 

The  computer-aided  analysis  was  repeated  using  the  Rexolite  rai.  : 
and  A  E-plane  patterns  were  computed  and  interpolated  as  before.  The 

htl-i 

suits  are  graphed  in  Figure  13  and  compared  to  the  measured  patterns  wit' 

the  actual  radome  in  place.  The  radome  was  oriented  with  respect  to  th>: 

antenna  such  that  the  radome  axis  of  symmetry  made  an  angle  of  1  “  c  with 

respect  to  the  antenna  z-axis,  and  the  tip  of  the  radome  was  located  in  ' 

$=-45°  plane  of  the  antenna  (Figure  2).  Agreement  is  good  for  the  primro 
E  :  I 

A  component  over  the  range  6  <36°;  the  discrepancies  outside  this  rr-s 
ELv  - 

of  angles  are  attributed  to  reflections  not  accounted  for  in  the  anaiys: 

E 

The  poor  agreement  between  the  cross-polarized  A  commorients  in  Facer'. 

ELp 

is  attributed  to  deficiencies  in  the  analysis,  particularly  the  absence  < 
any  accounting  for  reflections  inside  the  radome. 

DISCUSSION 

It  is  important  to  distinguish  the  band-limited  or  aperture-limit' 
nature  of  the  plane  wave  spectra  of  each  element  at  z=0  from  the  unlimit  i 
nature  of  the  PWS  of  the  four-horn  array  of  elements  at  z>0.  Each  ele rvi 
has  finite  size;  hence,  it  is  a  source  whose  aperture  fields  have  bounce 
support.  It  follows  that  the  PWS  is  not  limited  in  extent  in  the  wave  nut 
ber  plane  since  a  function  and  its  Fourier  transform  cannot  both  have-  fie 
support  [22]  .  The  unlimited  nature  of  the  element  plane  wave  spectra  rnu 
be  taken  into  consideration  when  implementing  the  extrapolation  algcnih: 
to  minimize  the  effects  of  aliasing. 


Tiic-  near  field  of  the  four-horn  array  can  be  synthesized  ire;  the 

element  near  field  or  from,  its  plane  wave-  spectra.  The  results  presented 

at  -ve  utilize,:  the  latter  method,  and  since-  no  windowing  function  was 

.i,  :  it.  the  plane  wave-  s;  ectra  of  the  four-horn  array,  t:.-.  respiting 

near  lit-Jds  are  tn--'-  -s-.  tied,  ly  those  rignt  at  the  z=c  plane  wner*  t:.v  *■  _  -  - 

me;,  t  s  are  locaiecj.  It  i  s  we  i .  known  tnat  trie  pdanar  tear  fie..; 

as  a  function  of  z[i  j;  viz.,  the  ilanc  wave  sj>octrar.  is  nodi  f :  ••  ;y  ‘ 

-  j  k  k  z  >.  ,  *■ 

exponential  function  e  °  “  sc-  that  E  (x ,  y , z)  =F-,  A  (k  ,k  )< 

n  f  x 

For  values  of  z  of  only  a  few  wavelengths,  the  effect  of  t:.e  e-xpcner.ti..  1 
function  is  to  drastically  attenuate  the  evanescent  mode s  of  the-  s:  •.  •*.  ra 
and  to  cause  the  planar  near  field  to  spread  out.  Thus,  what  was  a  near 
field  having  bounded  support  at  z=0  becomes  a  Fresnel  field  of  unlimited 
extent  with  increasing  z;  concurrently,  the  initial  spectra  having  un¬ 
limited  extent  become  limited.  In  the  limit  of  very  large  z,  the  spectra 
have  bounded  support  on  the  region  /kT^+  k^ t  1.  This  interchangeability 
of  the  roles  of  aperture-limited  PWS,  and  wavenumber-limited  near  fields, 
is  not  always  fully  appreciated. 

The  synthesis  procedure  described  above  also  allows  for  the  inclu¬ 
sion  of  virtually  any  amount  of  additional  measured  pattern  data.  The 
additional  data  form  added  constraints  on  the  near-field  solutions  finally 
obtained  in  the  same  way  that  the  principal  plane  data  have  worked  here. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  aperture  synthesis  procedure  described  has  been  demonstrated  to 
yield  excellent  results,  both  in  terms  of  the  measured  versus  synthesized 
patterns  and  in  the  reasonableness  of  the  near  fields  produced.  The  pro- 
is  specifically  applicable  to  computer-aided  radome  analysis. 


cedurc 


especially  in  regard  to  verification  of  cc-ru  sited  r«*«  ult  j  :  y  *r:  a:  . 
to  measured  pattern  data. 

It  is  recommended  that  additional  work  1-*  carried  out  \r  i  *. 
the  applicability  of  the  method  to  source  synthesis  using  com; lei*  -.r; 
tude  and  phase  data,  both  near-fielc  and  far-field.  Add i tier. a:  «our 
geometries  should  be  included  in  the  source  specification  to  ensure  c. 
conformance  with  physical  fact.  It  is  anticipated  that  the  actual  fi-.- 
of  a  number  of  common  source  Geometries  can  be  determined  quite  accura 
from  measured  radiation  data  using  the  general  procedure  described  her 
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INTRODUCTION 

This  paper  presents  the  salient  results  of  a  parametric  investiga¬ 
tion  of  radome  analysis  methods  [1].  Three  methods  of  analysis  were  in¬ 
vestigated  as  described  below.  Measured  boresight  error  data  obtained  at 
35  GHz  on  eight  combinations  of  three,  four-horn  monopulse  antennas  and 
five  polystyrene,  tangent  ogive  radomes  are  used  as  true  data  in  assessing 
the  speeds,  accuracies,  and  ranges  of  validity  of  the  three  methods  of 
analysis . 


METHODS  OF  ANALYSIS 


Three  computer-aided  methods  of  radome  analysis  were  investigated. 
Hie  theoretical  bases  for  all  three  methods  are  the  Huvgens-Fresne 1 
principle,  Lorentz  reciprocity,  and  geometrical  optics  [2].  The  first 
method,  called  herein  the  fast  receiving  method,  uses  geometrical  optics 
(ray  tracing)  to  fit  the  fields  incident  on  the  aperture  of  a  monopulse 
antenna  enclosed  by  ,_ne  radome  due  to  a  TEM  wave  incident  on  the  radome 
from  a  specified  direction  k  [3].  The  insertion  voltage  transmission  co¬ 
efficients  of  flat  panel  theory  for  parallel  and  perpendicular  polariza¬ 
tion  are  used  to  transform  the  plane  wave  fields  associated  with  each  ray 
from  their  values  at  the  incident  point  on  the  outside  surface  of  the  ra¬ 
dome  to  their  values  JE^,  H^  at  the  sample  point  in  the  aperture.  The 
voltage  received  by  each  channel  (F. ,  A  ,  A  )  t'ie  monoPulse  antenna 
is  obtained  according  to  the  reciprocity  integral 


'>  ‘  Js  da 

where  E  ,  H  are  the  aperture  fields  of  the  antenna  when 
where  S  is  the  aperture  surface. 


(1) 


transmitting,  and. 


The  second  method  of  analysis,  referred  to  herein  as  the  fast 
transmitting  method,  uses  a  transmitting  formulation  based  on  the  plane 
wave  spectrum  (TVS)  representation  of  the  antenna  aperture  fields  and  an 
equivalent  aperture  approach  [4].  The  aperture  fields  are  represented  by 
their  samples  at  FIN  equally  spaced  increments  in  x  and  y.  The  two-dimen¬ 
sional  Fourier  transform  of  the  aperture  fields  yield  the  corresponding 
piano  wave  spectra.  From  each  sample  point,  there  emanate  MN  plane  waves. 


Each  plane  wave  is  traced  to  an  incident  point  on  the  inner  surface  of 
the  radome  and  weighted  by  the  insertion  voltage  transmission  coefficients. 
The  modified  fields  of  each  plane  wave  are  then  added  together  at  each 
sample  point  to  produce  modified  aperture  fields  which  embody  the  effects 
of  the  radoiao .  The  voltages  received  by  the  radome-enclosed  antenna  are 
then  calculated  according  to  Equation  (1) ,  where  E^,  H  are  the  modified 
aperture  fields. 

The  third  method  of  analysis  is  referred  to  as  the  surface  integra¬ 
tion  method  [5].  The  voltage  received  by  the  antenna  enclosed  by  the  ra¬ 
dome  is  again  given  by  Equation  (1),  where  S  is  now  the  inside  surface  of 
the  radome,  _E-p,  H^.  are  the  radiated  fields  of  the  antenna  on  S,  and  E.,, 

Hp  are  the  fields  of  the  incident  plane  wave  on  S  as  transformed  from 
their  values  on  the  outside  surface  using  the  insertion  voltage  transmis¬ 
sion  coefficients.  For  some  directions  of  arrival  k,  portions  of  the  in¬ 
ner  radome  surface  are  "shadowed"  by  other  portions;  the  fields  on  the 
shadowed  portions  are  approximated  by  ray  tracing  as  described  above  for 
the  fast  receiving  case.  Aperture  integration  is  used  to  calculate  tae 
fields  of  the  antenna  at  each  point  on  the  radome  surface  from  the  soecied 
aperture  values. 

All  three  methods  were  implemented  in  Fortran  for  execution  on  tae 
Cyber  70  computing  system  at  Georgia  Tech.  Maximum  use  was  made  of  fea¬ 
tures  and  software  common  to  all  three  codes. 

ANTENNA/ RADOME  DESCRIPTIONS 

Three,  four-horn  monopulse  antennas  were  designed  and  fabricated 
exclusively  for  use  in  this  research,  corresponding  to  small  (BW-  ^=30°) , 
medium  (15°),  and  large  (8°)  sizes.  The  conical  horn  elements  of  each  an¬ 
tenna  were  machined  into  a  single  piece  of  aluminum.  The  salient  dimen¬ 
sions  of  the  antennas  are  given  in  Figure  1  and  Table  1.  More  details 
are  given  in  Reference  6.  The  elements  were  y-polarized. 

Five,  tangent  ogive  radomes  having  dimensions  given  in  Figure  1  and 
Table  2  were  machined  from  cylinders  of  polystyrene  (c  =2.54,  tan6=.OC2), 
corresponding  to  small,  medium,  and  large  radome  sizes.  A  fineness  ratio 
of  1.0  was  used  for  each  size;  in  addition,  for  the  medium  size,  fineness 
ratios  of  1.5  and  2.0  were  used.  All  five  radomes  had  a  wall  thickness 
of  0.25  inch,  corresponding  to  a  full  wavelength  in  the  dielectric  at  35 
GHz  and  a  design  angle  of  approximately  60  degrees. 

The  antenna  and  radome  were  used  together  in  eight  different  com¬ 
binations  for  measurement  and  analysis  purposes.  Hie  intent  here  was  to 
obtain  true  data  from  measurements  for  comparison  to  the  predicted  results 
over  ranges  of  parameters  which  would  help  clarify  any  deficiencies  in  the 
three  methods  of  analysis. 

RESULTS 

Comparisons  of  measured  and  computed  boresight  error  data  for  five 
an tenna/ radome  combinations  are  presented  in  Figures  2  through  5.  Bore- 
sight  errors  in  both  the  elevation  difference  channel  (BSEEL)  and  azimuth 
difference  channel  (BSEAZ)  are  shown  in  each  graph  as  functions  of  the 


radome  orientation  angle  a.  The  antenna  and  radome  were  mounted  together 
such  that  the  radome  axis  of  symmetry  passed  through  the  antenna  axis  of 
symmetry  at  the  gimbal  point  (Figure  lb)  at  an  angle  of  15  degrees.  The 
radome  could  then  be  rotated  about  the  antenna  axis  of  symmetry  through 
the  angle  a  indicated  in  Figure  la,  thereby  placing  the  tip  of  the  radome 
in  any  desired  plane  of  scan.  For  example,  a=0  corresponds  to  the  radone 
tip  being  in  the  azimuth  plane  (+xz-plane)  of  the  monopulse  antenna; 
ct=-90°  places  the  tip  in  the  elevation  plane  (-yz-plane  in  Figure  la) . 

In  Figures  2  through  5,  a  varies  from  zero  to  -90  degrees.  For 
ct=0 ,  symmetry’  dictates  that  there  should  be  no  boresight  error  in  eleva¬ 
tion,  and  the  azimuth  boresight  error  should  be  nonzero.  For  a=-90°, 
azimuth  boresight  error  should  be  zero,  and  elevation  boresight  error 
should  be  nonzero.  The  boresight  errors  are  defined  here  as  the  true  di¬ 
rection  to  the  target  in  the  antenna  coordinate  system  of  Figure  la;  e.g., 
positive  boresight  error  in  azimuth  (elevation)  means  that  the  target  lies 
in  the  +xz  (+yz)  plane. 

Figure  2  compares  the  computed  results  for  all  three  methods  with 
the  measured  boresight  error  data.  In  carrying  out  the  computations  for 
the  fast  receiving  case,  it  was  found  that  only  four  points  were  needed 
in  the  antenna  aperture,  located  at  the  center  of  each  element,  to  obtain 
essentially  the  same  results  as  were  computed  using  a  49  X  49  point  repre¬ 
sentation  of  the  aperture  fields  [7].  The  same  four-point  representation 
was  also  used  in  the  other  two  computer  codes  because  the  enormous  compu¬ 
tation  times  required  using  the  49  X  49  point  representation  were  prohi¬ 
bitive.  Even  using  the  four-point  representation,  the  surface  integra¬ 
tion  required  such  long  run  times  that  only  the  small  radome  could  be 
adequately  analyzed  on  the  Cyber  system. 

Examination  of  Figure  2  shows  that  none  of  the  three  methods  of 
analysis  accurately  predicts  the  measured  boresight  errors  in  the  case  of 
the  snail  antenna  with  small  radome,  especially  for  the  elevation  channel. 
'T'e  ‘nst  receiving  method  tends  to  overestimate  the  error  while  the  fast 
transmitting  method  underestimates  it.  The  surface  integration  method 
dot.  •  the  best  job  for  the  azimuth  error.  All  three  methods  predict  the 
correct  sign  of  the  error  for  most  of  the  range  of  u  considered.  Defi¬ 
ciencies  in  the  methods  of  analysis  to  account  for  the  lack  of  agreement 
w :  measured  results  cannot  be  isolated  on  the  basis  of  the  data. 

Figures  3  and  4  present  results  for  the  medium  antenna  and  medium 
-a domes  having  two  different  fineness  ratios.  The  fast  receiving  method 
predicts  the  measured  errors  most  accurately,  though  not  as  accurately  as 
desired.  The  fast  transmitting  method  underestimates  the  errors.  Un¬ 
stable  results  were  obtained  with  the  surface  integration  method  as  indi¬ 
cated  in  Figure  3,  indicating  unresolved  problems  with  the  computer  code. 
The  two  fast  methods  accurately  predict  the  sign  of  the  error,  and  the 
fast  receiving  method  does  predict  accurately  an  increase  in  errors  with 
-■creasing  fineness  ratio  as  expected. 


Figure  5  presents  the  results  for  the  medium  antenna  and  large  t.i- 
dome.  Neither  of  the  two  fast  methods  predicts  the  elevation  error  very 
accurately.  Better  predicted  results  are  obtained  in  the  azimuth  errot , 
v : th  the  fast  transmitting  being  more  accurate  than  the  fast  receiving 


The  computation  times  and  core  storage  recuirements  for  the  three 
methods  of  analysis  are  presented  in  Table  3  for  comparison  purposes.  tor 
the  fast  methods  of  analysis,  these  two  parameters  are  independent  of  the 
radome  size,  but  do  depend  on  the  number  of  points  in  the  antenna  aper¬ 
ture:  for  the  fast  receiving  method,  computation  time  increases  directly 
with  the  number  of  points;  for  the  fast  transmitting  method,  the  logarithm 
(base  2)  of  the  computation  time  increases  in  the  same  manner. 

The  computation  time  of  the  surface  integration  method  depends  on 
both  the  number  of  points  in  the  aperture  and  the  number  of  points  used 
to  represent  the  fields  on  the  radome  surface.  In  Figure  2,  four  points 
were  used  for  the  aperture  fields,  and  781  points  (spaced  A/3  apart)  were 
used  on  the  radome  surface.  In  Figure  3,  2291  points  at  A/3  spacing  were 
used  on  the  radome  surface,  and  four  points  were  used  for  the  aperture 
fields . 

Comparisons  of  measured  and  computed  values  of  on-axis  sum  channel 
gain  loss  relative  to  the  case  of  no  radome  are  shown  in  Table  4  for  three 
positions  of  the  radome  tip.  For  the  small  losses  encountered,  no  method 
consistently  predicted  the  measured  radome  losses. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  conclusion  to  be  made  from  the  data  presented  is  that 
there  is  room  for  improvement  in  the  predictive  accuracies  of  the  compu¬ 
ter  codes  used  with  respect  to  the  important  parameters  of  boresight  er¬ 
ror  and  radome  loss.  This  conclusion  is  especially  true  for  the  case  of 
small  antennas  and  radomes  where  the  effects  of  antenna /radome  interac¬ 
tions  are  not  properly  included.  For  the  case  of  moderate  sized  antennas 
and  radomes,  the  fast  receiving  method  is  attractive  because  of  the  fast 
computation  time  and  reasonably  accurate  results.  The  surface  integra¬ 
tion  code  used  requires  additional  development  before  any  valid  conclu¬ 
sions  can  be  drawn  about  its  predictive  accuracy;  however,  it  is  clear 
chat  its  applicability  will  be  restricted  to  small  radomes  because  of  the 
relatively  large  computation  times  required. 

It  is  recommended  that  the  fast  receiving  method  be  modified  so 
that  refractive  effects  and  ray  spreading  are  more  accurately  accounted 
for  in  the  ray  tracing  procedure.  It  is  also  recommended  that  the  sur¬ 
face  integration  method  be  further  developed,  even  to  include  first-order 
reflected  fields.  Additionally,  it  is  recommended  that  new  experimental 
techniques  be  developed  which  will  allow  the  determination  of  radome 
fields  close  to  the  dielectric  so  that  deficiencies  in  the  analytical 
methods  can  be  isolated  and  corrected. 
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Table  1.  Dimensions  in  Inches  of  Antennas.  (See  Figure  1(a)). 


An  tenna 
Identification 

Element 
Diameter  A 

Element  Kalf- 
Spacing  C 

Overall 
Diameter  B 

Small 

.1875 

.310 

.160 

.875 

Medium 

.250 

.620 

.320 

1.540 

Large 

.375 

1.225 

.878 

3.040 

Table  2.  Dimensions  in  Inches  of  Radomes.  (See  Figure  1(b)). 


Radome 

Identification 

Outside  Outside 

Diameter  D  Length  C 

Inside 
Diameter  E 

Inside 
Length  B 

Small  (F=l .0) 

2.55 

2.35 

2.05 

2.05 

.437; 

Medium  (F=1.0) 

4.00 

3.81 

3.50 

3.50 

.75 

Medium  (F=1.5) 

4.00 

5.66 

3.50 

5.25 

.75 

Medium  (F=2.0) 

4.00 

7.52 

3.50 

7.00 

.75 

Large  (F=1.0) 

6.91 

6.72 

6.41 

6.41 

.375 

Table  3.  Computation  Time  and  Core  Storage  Requirements. 

Computation  Time *( se c. ) /Core  Storage  (octal) 


Case 


Fast 

Receiving 


Fast  Surface 

Transmitting  Integration 


Small  Antenna,  Small  Radome 
Medium  Antenna,  Medium  (F=1.0) 


.24/67100  .62/77300  34.0/135300 

.25/67100  .62/77300  44.0/135300 


*Per  look  angle. 


Table  4.  Radome  Loss  Comparisons. 


Antenna 

Radome 

Alpha 

(Deg.) 

Measured 

Loss 

Fast 

Revg 

(dB) 

Fast 

Xrct-E 

Surface 

Integration 

Small 

Small 

0 

0.4 

0.5 

0.5 

1.9 

M 

M 

-45 

0.7 

0.5 

0.4 

1.7 

n 

I! 

-90 

0.6 

0.3 

0.3 

1.9 

Medium 

Medium 

(F=1.0) 

0 

1.4 

0.4 

0.6 

1.2 

11 

*» 

It 

-45 

1.4 

0.4 

0.5 

1.2 

1 1 

M 

11 

-90 

0.8 

0.3 

0.4 

1.2 

Medium 

Medium 

(F=2 .0) 

0 

0.2 

0.1 

0.6 

- 

11 

f  r 

II 

-45 

0.5 

0.1 

0.4 

- 

♦  1 


-90 


0.6 


0.1 


0.2 


(a)  Antenna  Geometry  Showing  Radome  Orientation 


B 


C 


I 

h)  Hulum*  Geometry 

f  icjuri*  1  Antenna  and  Raclome  (ieom»-tn**s  .»•'  Dinien- 
■n  Tables  1  arul  2. 
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